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Capítulo 1 
Resumen 



Éste documento tiene como objetivo hacer una introducción a la teoría 
de los monopolos magnéticos clásicos y cuánticos, introduciendo la electro- 
dinámica desde el punto de vista conceptual intentando hasta donde sea 
posible no descuidar la formulación matemática con el fin de explicar tanto 
física como matemáticamente las ecuaciones de Maxwell; posteriormente se 
hará una breve pero profunda exposición de las bases de la relatividad espe- 
cial que tiene como objetivo usarla para que junto con la electrodinámica se 
pueda formular la teoría covariante de la electrodinámica y analizar su inva- 
rianza de Lorentz. El documento busca modelar el campo electromagnético 
del monopolo magnético con el objetivo de compararlo con el campo elec- 
tromagnético de los polos eléctricos estáticos, y en movimiento. A manera 
introductoria se puede observar que el campo de los monopolos es interesante 
pues en la electrodinámica clásica se ha supuesto la no existencia de los po- 
los magnéticos libres (PM). Clásicamente los monopolos se definen mediante 
una rotación de 7r/2, lo que se llama una transformación de dualidad. Dirac 
[1], demostró que el polo magnético sería necesario para cuantizar la carga 
eléctrica, pero se encontró un problema muy serio con una singularidad para 
6 = (O, Tc) que se llamó la cuerda de Dirac. Aunque nunca hubo evidencia 
experimental del monopolo y Dirac al final de su vida puso en duda la exis- 
tencia de éste, parece ser que la falla está en elementos erróneos de la teoría 
de Dirac. Adicionalmentelos monopolos relacionan muy cercanamente la físi- 
ca de partículas y la cosmología, pues a muy altas energías se restablece el 
grupo de unificación, con cuyo enfriamiento la simetría fue espontáneamente 
rota. También la teoría puede definir un nuevo tipo de interacciones a un or- 
den tal que se podría cambiar el número bariónico llevando a la predicción de 
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CAPÍTULO 1. RESUMEN 



que el protón sería inestable. Finalmente en la teoría de unificación habrían 
partículas estables que llevan cargas magnéticas o monopolos magnéticos. 
Una de las teorías postuladas para la unificación es la teoría de las cuerdas, 
aunque ésta aún está en desarrollo es muy matemática y en este momento, 
no parece que sea posible su verificación experimental lo que la hace una 
teoría que sea de mucho agrado para los fíeos. 

Las críticas a los monopolos son muy interesantes [2], [3], [1], [5] pero 
solo desde el punto de vista de la simetría de las ecuaciones de Maxwell, es 
innegable la importancia del monopolo magnético. La dificultad para estudiar 
experimentalmente al PM es que se requerirían aceleradores con una energía 
del orden de 10^^ GeV, los cuales aún no se han construido. 



Capítulo 2 

Introducción a la 
Electrodinámica 



Las ecuaciones de Maxwell en unidades Gaussianas fiQ = Eq = 1 (para 
convertirlas al sistema MKS ó SI, se debe tener en cuenta la equivalencia 
£o = y ¿"o = 47r X 10-^)son: 

IdD 4n 
V.D = 47rp Vxíf — = —J 

^ c dt c 

1 dB 

V.B = Q V X E + -—^0 (2.1) 

c ot 

en las cuales la Ley de Gauss, V . Z) = V . -E = Anp, significa que la 
densidad de flujo eléctrico diverge desde un foco (densidad de carga eléctrica 
positiva) o hacia un sumidero (densidad de carga eléctrica negativa), con lo 
que se pueden graflcar las líneas de campo eléctrico de forma radial. Así en 
una región del espacio tridimensional se construye una superficie cerrada 
imaginaria, denominada superficie Gaussiana, con el objetivo de analizar el 
flujo de las líneas de campo. Si se considera el espacio en un volumen r, la 
Ley de Gauss toma la forma, 

J^iV .D)dT^ÍTT J^pdr, (2.2) 

se puede observar que hay una relación directa entre ese volumen y una 
superflcie cerrada por medio del Teorema de la Divergencia o de Gauss, 
quedando 

D.dS^ inq. (2.3) 
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CAPÍTULO 2. INTRODUCCIÓN A LA ELECTRODINÁMICA 



Si el campo es constante, para una carga puntual, que es simétrica con una 
superficie esférica se llega a la conocida Ley de Coulomb, 



que el sistema MKS es 



D-^,r (2.4) 



^ = ir-^- (2-5) 



donde se usó la relación D = SqE, que relaciona los dos campos vectoriales 
eléctricos (el campo densidad de flujo eléctrico D j e\ campo intensidad de 
campo eléctrico E). 

Teniendo en cuenta la identidad del cálculo vectorial Vr~^ = — r~^r, se 
puede reescribir la Ley de Coulomb quedando 



E^-V [ — ^ = -V* (2.6) 

donde $ es una función escalar, o campo escalar denominada potencial 
eléctrico, que define el trabajo por unidad de carga para mover ésta de un 
punto a otro en contra del campo eléctrico. 

La Ley de Amp ere- Maxwell, V x ií" — \^ — dice que si existe 

una densidad de corriente eléctrica, habrá una circulación de la intensidad 
de campo magnético alrededor de ésta, pero para poder conservar la conti- 
nuidad de corriente eléctrica, es necesario que la densidad de fiujo eléctrico 
sea temporalmente variante. Lo anterior se puede observar analizando que la 
divergencia de un rotacional siempre es cero, piics el rotacional produce una 
circulación, que no posee divergencia. Por tanto si existe variación temporal 
del campo eléctrico existirá un campo magnético sin necesidad alguna de 
fuentes. Así, 

lo que se reduce a 



c \ dt c 



- ^ V . D = 47rV . J (2.8) 

dt ^ ^ 

y mediante el uso de la Ley de Gauss, finalmente se llega a la ecuación de 
continuidad. 
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que no se verificaría si no hubiera variación de la densidad de flujo eléctrico. 
Adicionalmente debido a la circulación del campo magnético, éste está defi- 
niendo una trayectoria cerrada o una superficie abierta. Si la Ley de Ampere- 
Maxwell se integra sobre una superficie abierta, se llega a 



1 d í An f 

(V xH).dS--— / D.dS^— / J.dS (2.10) 

^ ^ cdtJs c Js ^ ^ 



en donde se puede observar que Jg J . dS es una corriente eléctrica, que es 
originada por un movimiento de cargas puntuales en una dirección determi- 
nada y local o sea mediante un vector densidad de corriente eléctrica. Hay 
que hacer claridad que en éste punto se definen las tres corrientes eléctricas 
de la electrodinámica; la corriente de conducción que se origina en la Ley 
de Ohm, J = crE, donde a es la conductividad; la corriente de convección, 
J = pvd, donde p es la densidad volumétrica de carga y es la velocidad de 
deriva o de arrastre que depende del tipo de conductor por el que se mueven 
las cargas libres, y la corriente de desplazamiento de Maxwell, Ja = 
que se origina con la variación temporal de la densidad de flujo eléctrico. 
Usando el Teorema de Stokes, la Ley de Ampere- Maxwell integral es 

H.dL^—I (2.11) 

c c 

donde la corriente / contiene tanto la Ley de Ohm como la corriente de 
desplazamiento de Maxwell. 

La ley de Faraday-Henry, Vx£^-|-i^ = 0, informa que si la densi- 
dad de flujo magnético es temporalmente variante, se originará una contra- 
circulación de la intensidad de campo eléctrico alrededor de ésta. Esta es la 
base de la Ley de Lenz y del origen de las corrientes parásitas, o corrientes 
de Eddy. Por tanto una variación temporal del campo eléctrico produce una 
ciculación levógira del campo magnético alrededor de ésta. Realizando una 
integración sobre una superflcie abierta se llega 

/ (V xE).dS^--^ í B.dS^--^ (2.12) 
Js cdtJs c dt ^ ^ 

donde utilizando el teorema de Stokes se veriflca /^(V x E) . dS = §^E) . 
dL — fem y por tanto, 

1 9$ 

f8m = ---^. (2.13) 
c OI 
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CAPÍTULO 2. INTRODUCCIÓN A LA ELECTRODINÁMICA 



donde fem es la fuerza electromotriz, que no tiene nada que ver con una 
fuerza pues tiene unidades de voltios, pero se sigue llamando así, aunque no 
sea correcto, debido a 'razones históricas'. Es de observar que la derivada 
temporal del fujo magnético en la Ley de Faraday-Henry (ec. (12.121) ) afecta 
tanto al campo como a la sección areal; por tanto. 



donde en el primer término del segundo miembro la sección areal es cons- 
tante y el campo varía temporalmente y en el segundo término el campo es 
constante y la sección areal varía respecto al tiempo. 



Finalmente, la llamada ley de Gauss del magnetismo, V . S = V . ií = 
O pues /io = 1, muestra que el campo magnético no es divergente y por tanto 
los polos magnéticos (fuentes magnéticas) libres no existen. Si se analiza en 
todo el espacio, con un volumen r. 



usando el teorema de la divergencia, dicha ley se escribe en su forma integral 
como 




(2.14) 



Problema: Describa un sistema físico en el cual estén 
presentes los términos del segundo miembro de la Ley 
de Faraday-Henry. 




(2.15) 




(2.16) 



2.0.1. Potenciales Vectoriales 



Si se consideran detenidamente las ecuaciones de Maxwell, 



1)V .D 



4vrp, 




3) V.fí 



O 



4:)V X E 



IdB 



O 



(2.17) 



+ - 



c dt 
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se puede observar en la ley de Gauss del magnetismo (ausencia de monopolos 
magnéticos libres, Ec. 3) de (12.171) ) el hecho de que la divergencia del campo 
magnético es nula, entonces dicho campo está definido por el rotacional de 
otro campo al que se da el nombre de potencial magnético, 

V .B = V .V X A = (2.18) 

ésto es B = 'V X A. Usando de la Ley de Biot-Savart, 



, 1 f J(X') X (X — X') ,n , 

B(x) = - / -^—^ ^— ^áV 2.19 

^ ^ c x-x'^ ^ ' 



se nota claramente una analogía con la Ley de Coulomb, 



N / p{x')ix — X') ,o , , , 

E{x) = I ^ dV (2.20) 



\X — X' 



lo que lleva a dilucidar la forma matemática del potencial vectorial. Éste es 

1 í Jix' 



Aix) = - / , ^ \ ciV (2.21) 
^ ' c J \x - x'\ ^ ' 

pero dicha definición llevaría a ajustar al potencial en un punto específico del 
espacio y no en cualquier punto. Por tanto a la ecuación (12.211) se le adiciona 
el gradiente de una función escalar V^(íc) de forma tal que se calibre al 
potencial vectorial de la forma A A + V^, lo que se llama realizar una 
transformación de gauge en la teoría, y lo que se busca con ésto es demostrar 
que la teoría es invariante de gauge. También para el potencial escalar se 
verifica una transformación de gauge análoga, $ — > $ — Si se reemplaza 
S = V X A y, adicionalmente se usa la ec. (12. 6p . en la ecuación 2) de (12.171) . 
se obtiene, 

V X (V X A) - 7 = —J (2.22) 

c ot c 

donde usando identidades vectoriales se puede modificar a 

V(V.A)-V^A + Í^ = ÍÍJ (2.23) 

c ot c 

pero se observa que hay algo que no funciona, pues si hay campos variantes 
con el tiempo, la ecuación 4) de (I2.17P queda 

„ _ 19(VxA) „ / ldA\ 
VxE + — ^— ^ = VxLE + -— =0 (2.24) 



c dt \ c dt I 
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CAPÍTULO 2. INTRODUCCIÓN A LA ELECTRODINÁMICA 



por lo que el campo eléctrico queda definido no como la ec. ( 12.6p . sino como 
E + ^^ = -V$, y por tanto E = -V$-i^. Realizando éstos reemplazos 
en la Ley de Faraday-Henry, se obtiene una ecuación que posee la forma de 
una ecuación de onda, pero con potenciales acoplados: 

C Ot C 

V(V.A)-V^A + -V- + -^ = -J 

Para desacoplarlos se creó la condición de Lorenz, 

1 9$ 

V.A+-— = 0, (2.26) 

c ot 

quedando la ecuación de onda para el potencial vectorial en presencia de 
fuentes de campo magnético, 

Usando las transformaciones de gauge de los potenciales en la condición de 
Lorenz, 

V.(A + V*) + -5^(*-- 

se observa que la función escalar \E' con la que se calibra la electrodinámica, 
debe cumplir la ecuación de onda 

1 d'^'^ 
c ot^ 
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Finalmente reemplazando el valor de la intensidad de campo eléctrico E 
_V$ - en la Ley de Gauss V .E = Anp, 

c~dt 



V . 



-Vi 



Airp 



1 dA 

V . V$ + V . — — = -Ano 
c ot 



1 



(2.30) 



donde se usó la condición de Lorentz (ec. 12.261) . El resultado final es que los 
potenciales escalar y vectorial se comportan de forma ondulatoria y que las 
ecuaciones de Maxwell se pueden reescribir en función de éstos, y no de los 
campos eléctrico y magnético. 

Se observa adicionalmente que los campos E y B sin presencia de fuentes en 
el vacío satisfacen la ecuación de onda, 



V^E = 



1 d^B 



[2.31) 



y por tanto usando la teoría de Fourier, el campo eléctrico toma la forma 
1 



Eir.t) 



i2n 



Ei{k)e 



ik.r—iujt 



+ E2Í-k)e 



ik.r+iuit 



d-'k 



(2.32) 



y el campo magnético toma una forma similar, 
1 



B{r,t) 



i2n 



Bi{k)e 



ik.r—iuit 



+ B2Í-k)e 



ih.r+iu)t 



(2.33) 



donde ^ = k = \k\. siendo k el número de onda, y donde el carácter ondula- 
torio de los campos electromagnéticos se hace evidente. 



2.0.2. Simetría de las Ecuaciones de Maxwell 

A pesar de que éstas ecuaciones (ver ec. ( 12. ip ) nos dan la idea de que la 
electrodinámica no es simétrica en presencia de cargas, se puede notar que 
desde el punto de vista energético si existe tal simetría, pues la densidad de 
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CAPÍTULO 2. INTRODUCCIÓN A LA ELECTRODINÁMICA 



energía electromagnética es u = {8'7i)~^{E.D + B.H), donde es evidente que 
las densidades de energía eléctrica y magnética son simétricas entre si, lo cual 
es aún más notorio en la forma matemática que tiene el flujo de radiación 
S = c{Att)-\E X B). 

Si se supone la existencia de fuentes magnéticas, las ecuaciones de Ma- 
xwell (ver ec. 12.11) se reescriben de la siguiente forma: 

V .D = 47rpe V X H 



IdD 




c dt 


C 


IdB 




c dt 


= — Jh 

c 



V.B = Anph -V xE--— = —Jh (2.34) 

c ot c 

en las cuales es evidente tal simetría, donde se puede observar que aparecen 
las fuentes magnéticas (polo magnético g con su densidad de carga nagnética 
asociada ph y la densidad de corriente magnética Jh)- Aunque la motivación 
primordial que se expone aquí es la simetría, la existencia del monopolo es 
requerida en las teorías que estudian la estructura de la materia a escalas 
muy pequeñas (del orden de 10~^° m) y para poder experimentar allí, las 
energías necesarias son tan grandes que no existen aceleradores para hacerlo, 
lo cual llevó a algunos científlcos a trabajar en cosmología, pues en ésta área 
se encuentran cuerpos y procesos en la naturaleza con mucha energía, tal 
como explosiones estelares, pulsares, etc., mientras otros buscan el desarrollo 
teórico de los monopolos, desde el punto de vista de la uniflcación, teorías de 
cuerda, o teorías supersimétricas. 

Inicialmente Dirac creó el monopolo magnético para cuantizar la carga eléctri- 
ca, pero actualmente se aplica a las interacciones fuertes y débiles, a las 
teorías de gauge abehanas y no-abelianas, también para el estudio de los 
procesos de conflnamiento de quarks en la cromodinámica cuántica (QCD). 
Actualmente el problema tiene que ver con el estudio de la singularidad 
intrínseca que tiene la teoría de monopolos de Dirac, o cuerda de Dirac, 
que al parecer debido a que nunca se ha medido y a que todos los teore- 
mas matemáticos fallan pues sobre la superflcie abierta se crea una línea de 
singularidad, se concluye que hay una falla conceptual en la teoría de Dirac. 



Capítulo 3 
Relatividad Especial 



Debido a que se pensaba en los años 1900 que la luz necesitaba un me- 
dio para propagarse, pues el sonido necesitaba un medio para hacerlo, las 
olas necesitan para propagar su perturbación un medio fluido, etc., los físicos 
propusieron una substancia que llenaba el espacio, que permitía que la luz 
se propagara en el vacío, y lo llamaron ctcr. Las propiedades de éste ctcr 
son: no tiene densidad, no tiene alguna propiedad fisico-química u óptica, 
pero produce un viento cuando la tierra se mueve a través de él, lo que per- 
mite determinarlo midiendo la velocidad de la luz en el marco de la tierra 
en posiciones tales que se puedan cumplir las transformaciones de Galileo 
para que al recombinar los haces de luz se observen procesos de interferen- 
cia. El dispositivo para realizar éste experimento fue llamado interferómetro 
de Michelson-Morley. Como se mencionó anteriormente se basaron en las 
transformaciones de Galileo que correlacionaban dos sistemas de coordena- 
das inerciales bajo la medición de un evento en el espacio; éstas son, 

r'^r-vt, t'^t. (3.1) 

Debido a que los resultados del experimento de Michelson-Morley fueron 
negativos, aiin con luz extratcrrcstre, el resultado no fue completamente en- 
tendido hasta que apareció la teorí a de la relatividad especial de Albert 
Einstein en la cual propuso dos postulados: 

Postulado de la relatividad. En éste postulado Einstein dijo que las 
leyes de la naturaleza y también los experimentos realizados en un marco 
de referencia incrcial específico son invariantes con respecto a cualquier otro 
marco de referencia inercial. 
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CAPÍTULO 3. RELATIVIDAD ESPECIAL 



Constancia de la velocidad de la luz. La velocidad de la luz es inde- 
pendiente del movimiento de la fuente que lo origina. 

Einstein se enfrentó a un problema cuando al observar las transforma- 
ciones de Galileo (ver ec. 13.11) . se dió cuenta que la ecuación temporal era 
sospechosa, pues pensaba que el tiempo no era invariante entre los sistemas 
de coordenadas primado y no-primado, por lo que se debían cambiar en lugar 
de las ecuaciones de Maxwell (ver ec. 12.11) que eran existosas experimental- 
mente bajo todo punto de vista. 

Para aplicar éstos postulados, se crea un experimento en el que existen dos 
sistemas de referencia inerciales primado (SRP) y no primado (SRNP). El 
SRNP está estático y el SRP se mueve con velocidad constante, paralelo al 
eje X en el sentido positivo. Cuando coinciden los dos orígenes, en el SRNP 
se enciende un foco y la luz se propaga en tres dimensiones a través del es- 
pacio. Los dos sistemas verifican las ecuaciones (si se desea profundizar más, 
consulte los textos de Jackson [6], Landau [7] y Goldstein |9]), 



+ + z"^ - (ct)^ = O 
{x'f + {y'f + {z'f - {cff 



O 



(3.2) 



en donde se usó el segundo postulado en el cual c = c'. Debido a la iso- 
tropía del espacio y al primer postulado se concluye que las ecuaciones de los 
dos sistemas deben ser iguales, y que la diferencia entre ellas incide en una 
constante multiplicativa que marca el cambio de escala entre los marcos de 
referencia. Por tanto se tiene que. 





lx'\ 




/Aii 


A12 


Al3 


Al4^ 




íx\ 




y' 




A21 


A22 


A23 


A 24 




y 




z' 




A31 


A32 


A33 


A 34 




Z 




V/ 




\A41 


A42 


A43 


A44y 




vJ 


debido a la isoptropía del espacio A 


12 = 


Al3 = 


= A21 




A23 


A34 — A42 — A43 — 0, 


y por tanto. 










íx'\ 












Al4\ 




íx\ 




y' 







A22 










y 




z' 










A33 







z 


Vi 












A44/ 




VV 



^24 



A 



31 



\32 



Las ecuaciones deben ser lineales y en v = O deben coincidir. Por tanto An = 
A22 = A33 = A44 = 1, y ahora si se considera el movimiento paralelo al eje x 
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teniendo en cuenta que debe existir una similitud con las transformaciones 
de Galileo y por tanto A14 = —vXu (ver ec. I3.ip se llega a 



/x'\ Aii O O -vXu\ fx\ 

y' _ O 1 O O y 

z' ~ O O 1 O z 

\t'J \A4i Xu J \tj 



Reemplazando éstas ecuaciones en la ecuación principal x"^ + y'^ + z'^ — (cí)^ = 
{x'f + {y'f + {z'Y - {ct'f se obtiene 

x^ + y'^ + z^ - [ctf = [Xu{x - vt)]"^ + y'^ + z^ - [c{Xiix + Xi4)f . (3.3) 

Al resolver éste sistema de ecuaciones, y con /3 = -, se llega a las transfor- 
maciones de Lorentz (TL), 



X 



t' 



X — vt 



y =y 



z = z 



t-P-x 

ti 



(3.4) 



Problema: Demuestre las ecuaciones (13.41) mediante el 
uso de la ecuación (13.31) . 

Usando las notaciones /3 = ^, 7 = (1 — /9^)~^/^ y xq = cí, xi = x, X2 = 
y, X3 = z, se reescriben las TL como 

Xo = 1 (xo - pxi) 
x[ = 7 (xi - Pxq) 

X2 — ^2 

X3 = xs. (3.5) 



Si el SRP se mueve con una velocidad en 71^, se puede ver que las ecua- 
ciones se escriben de forma diferente, pues la cantidad /5 debe ser un vector y 
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por tanto los productos comunes entre escalares se convierten en productos 
punto. Para la parte espacial se tiene que, 



X 1 



Xn 



1 + 



7 



Pt X2 + 



P' 

7-1 



1 + 



7 



Pt X, + 



p 

7-1 



ÍP1P2X2 + P1P3X3) - iPlXo 

2 {P2P1X1 + P2P3X3) - 7P2X0 



{P3P1X1 + P3P2X2) - 1P3X0 



p2 ^-ly - p2 

y por tanto al reunir éstas ecuaciones se construye la ecuación vectorial 

7-1 



(3.6) 



X = X 



P' 



{f3 .x) f3 - -ff3xo. 



Para la parte temporal, 



Xn 



7 {^0 - Pi^i - P2X2 - P3X3) 



que vectorialmente se reescribe, resultando 



Xn 



j {xq - f3 . x) . 



(3.7) 



(3.8) 



(3.9) 



Problema: Demuestre las ecuaciones fl3.7p y (13. 9p me- 
diante el uso de la ecuación (13.51) . 



Teniendo en cuenta que la relatividad está definida en el espacio-tiempo, 
se puede escribir un cuadri-vector de la forma x'^ = x'^^lx^, x^, x^, x^), para fi 
0, 1, 2, 3, y usando la métrica g'^'^ = g^^ tal que 



9'" = 9,. 



/l O O O \ 
0-10 o 

0-10 

\0 o 0-1/ 



(3.10) 



donde se verifica que g^ug^^ = 9^ = hxi, Y en forma general que g^a9°"^ = 



d,A^ = ^ + V. A, 



Usando las notaciones = 9^ = d^Afj_ 

Q 02 — V , se construye una transformación relativista entre dos 
sistemas de coordenadas de la forma, x'^ = K^^^x^ + a'', (ésta que es análoga 
matemáticamente a la transformación de Galileo en la cual r' = r — vt, 
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solo que en la transformación relativista hay una rotación espacio-temporal, 
así como una traslación espacio-temporal) la que expresa propiedades que 
le dan el nombre de Grupo de Poincaré (GP). En el caso de la relatividad 
especial se toma el caso en el que a'^ = 0. La anterior restricción convierte el 
GP en el Grupo de Lorentz (GL) que son dos grupos que obedecen el álgebra 
de Lie, y en donde se usa la condición. 



Adicionalmente se verifica que 



X 



(3.11) 



(3.12) 



lo que significa que dos transformaciones sucesivas son también transforma- 
ciones de Lorentz y por tanto el álgebra es cerrada. La matriz de trans- 
formación A, o boost de Lorentz, se puede construir usando los postulados 
fundamentales de la relatividad especial y luego aplicándolos a un sistema 
tridimensional. Usando ¡3' = el boost de Lorentz queda definido por la 
ecuaciones ya anteriormente demostradas (ver ec. (13.71) y (13. 9p ): 



r' = Ar 

donde los vectores r', y r, son vectores columna definidos por 



(3.13) 



x[ 



y por r 



Xi 
\X3/ 



(3.14) 



y donde el boost de Lorentz, A, se define como 



A 



V- 



7 


-7/3i 


-7/^2 


7/5i 


1 + 


/3'/3i/32 


7/52 


/?'/?l/32 


1 + (3' (31 


7/?3 




/3'/32/?3 



-7/^3 

/3'/3i/33 
1 + P'PIJ 



(3.15) 



Problema: Demuestre la ecuación (I3.15p . 



El tiempo propio de cada sistema, que marca su evolución espacio-temporal 
y que es el tiempo tal cual se ve en el marco inercial del sistema, está definido 
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CAPÍTULO 3. RELATIVIDAD ESPECIAL 



por dr"^ = c^dt^ — dx^ = —g^^dx'^dx" , donde se constata que en otro sistema 
de referencia primado, el tiempo propio definido por 

dr''^ = -g^ydx^dx" = ~9^iu-g^-Q^dx"dx'^ 

= dr^ (3.16) 

es invariante bajo transformaciones de Lorentz. En la anterior ecuación se 
demuestra adicionalmente que gaf3 = g^u^^^^- Si se aislan las partes es- 
paciales y temporales: dx'^ = A\cdt y cdt' = A^^cdi y se dividen entre sí, se 
llega a 

entonces, 

A'o = AA% (3.18) 

lo que muestra una relación entre la parte espacio-temporal y la parte neta- 
mente temporal del GL. 

Si en la ecuación (13.110 . se imponen los valores a = /? = O entonces 
A^^A^Qg^^ = goQ = —1 donde separando la parte espacial y la temporal 

3 

la ecuación se transforma a " (A°o)^ = ~1 Y usando el resultado 

i=l 

hallado anteriormente: A*o = PíA^o^ tiene (/?^ — 1)(A°q)^ = —1, que al 
despejar se obtiene, 

A% = J—^ = 7 y A\ = pa- 

También se observa que A^-A^-g^^ = 1, donde siguiendo un procedimiento 
análogo al anterior resultan las ecuaciones, 

A°,=7/5, (3.19) 

que sirven para verificar el boost de Lorentz (ver ec. (13.151) ). Finalmente se 
escriben los vectores, 

x'o = lixo -13.x) 

x' = x + (/3 . x)l3 - 7/3x0, (3.20) 
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que marcan la estructura general de las transformaciones de Lorentz entre 
dos sistemas de referencia que se mueven en cualquier dirección. 

Problema: Realice un análisis físico de las ecuaciones 
(13.201) y halle las condiciones físicas para las cuales éstas 
son relevantes tanto para la velocidad como para la ace- 
leración. 



Con éste resumen de la relatividad especial y del grupo de Lorentz, se se 
sientan las bases para poder calcular la invarianza de la carga eléctrica y la 
covarianza de la electrodinámica en la siguiente sección. 



CAPÍTULO 3. RELATIVIDAD ESPECIAL 



Capítulo 4 



Invarianza de la Carga 
Eléctrica y Covarianza de la 
Electrodinámica 



Esta sección es netamente educacional, y está extraida primordialmente 
del texto de J. D. Jackson ¡6] y apoyada en los textos de L. D. Landau [7j y 
de H. Goldstein |9]. 

La importancia de hallar ésta invarianza es que ésta demuestra que la 
electrodinámica y las fuentes electromagnéticas respetan la relatividad espe- 
cial. 

Teniendo en cuenta la fuerza de Lorentz, aplicada a una partícula de carga 



donde el cuadri-momento está definido por = {pq = -f,p) = ""^(t^o = 
7c, U = ■yv), y donde se observa que la parte temporal es proporcional a una 
energía y transforma según la ec. fl3.18p . 
Adicionalmente, 




(4.1) 




(4.2) 
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y por tanto, = . E = j^, que en notación tensorial usando E .U = 
F^^U^, queda como ^ = |F°'^?7^. Considerando la densidad de corrien- 
te como J{x,t) = ^„e„5^(a3 — a;„(í))^^p^, y la densidad de carga como 
J° = en{S^{x — Xn)), donde 6{x) es la función de distribución delta de Dirac 
(ver el apéndice 19.11) . construimos la cuadri-densidad de corriente J"(x) = 
^^e„5^(íc — x„(í))^^^^, y por tanto la cuadri-divergencia estará dada por 



V . J{x,t) 



dxjjt) 
dt 



d 



- ^ e„— - Xn{t)) 



'dt 

n 

de{x, t) 



(4.3) 



que en cuatro dimensiones queda como = 9q,J" = O y de ésta manera la 
cuadri-corriente es invariante de Lorentz. Ahora definiendo la carga como la 
parte temporal de la cuadri-densidad de corriente, Q = J (PxJ^{x) se puede 
analizarla en función de todo el espacio como una variación temporal de la 
forma ^ = / d^x ^"^q^o^ ■ Usando el teorema de la divergencia se transforma la 
integral volumétrica a una integral sobre la superficie que encierra el volumen: 
^ = / d^xV . J{x) = O, la que se anula. En cuatro dimensiones ésta carga 
queda como Q = J d^xJ'^{x)da6{i]¡3X^) donde 6{x) es la función escalón o de 
Heaviside definida por 




X > O 
X < O 



y ?7/3 se define como r/i = 772 = ^73 = O, tjq = +1. Realizando una transfor- 
mación de Lorentz sobre Q resulta Q' = J d'^xJ'^{x)daO{r]'i^x^) y por tanto 
= A°'p'r]a. Si se calcula la diferencia entre las cargas transformadas me- 
diante la TL, se obtiene (Q' - Q) = {Q' - Q){e{rj'^x^) - e{r]px^)) = O pues 
la cuadri-densidad de corriente J" tiende a cero cuando |x| tiende a infi- 
nito y 6{rj'pX^) — Oijjpx^) tiende a cero cuando el tiempo tiende a infinito 
con X fijo. Por tanto Q es un escalar y es entonces invariante bajo cualquier 
transformación. 



25 



Las ecuaciones de onda para los potenciales A y están dadas por, 



1 d^A 

1 



AttJ 

c 

—Airp 



restringidas por la condición de Lorentz, 

-— + V.A = 0. 

c ot 

y cuyas soluciones tienen la forma de Fourier, 



(4.4) 



(4.5) 



f{r,t) 



dujfoj{r)e 



-iujt 



1 

/u,(r,í) = — y dtf{r,t)e^' 



(4.6) 



Si se define el cuadri-potencial A"' — ($, A) = (A^, A^), (para ¿ = 1, 2, 3), las 
ecuaciones de onda se mezclan dando como resultado, 



c 



0. 



(4.7) 



(4. 



En la electrodinámica los campos se definen mediante los potenciales de 
la manera usual, 

V$ 

c dt 

B = VxA 

que en forma de componentes quedan, 

Id Ai 9$ _ 

dxi 

= -{d'^A' 



c dxo 
dAk 



-(9^1 - 9M") 

d^A'') 



(4.9) 



donde se tuvo en cuenta que d"' — (g|^, —V). Con la anterior notación se 
escribe el tensor electromagnético F^^ — d'^A^ — d'^A'^, 



/o 


-E, 


— £'2 









-Bs 


B2 


E2 


Bs 





-Si 


\E3 


-B2 


Bi 


J 



(4.10) 
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donde el tensor electromagnético totalmente covariante está descrito como 



í 





El 


E2 


Es\ 




-El 





-B3 


B2 






B3 





-Bi 


\- 


-E3 


-B2 


Bi 


/ 



(4.11) 



Problema: Demuestre las ecuaciones fl4.1üp y fl4.1ip . 



Así, queda definida la transformación entre tensores F'^'^ — > F^^, defi- 
nidos en el espacio normal. En el espacio dual también se define el tensor 
electromagnético dual 



■fiu 



7¿ 



/O -Bi -B2 -B,\ 

Bi O £-3 — E2 

B2 — E^ O El 

\53 E2 -El o / 



donde el tensor e^'^'^^ está definido como, 



' + 1 para (a,/?, 7, 5) = (0, 1,2,3) 
permutación par 
— 1 para cualquier perm. impar 
^ O cualquiera de índices son iguales 



con la propiedad e^J,u'^l, = 
Dos de las Ecuaciones de Maxwell (Ec. 12.11) son 

dEi 



dxi 



4tt 47r 
47rp = — cp = — Jo 
c c 



_l9^_4vr 
C 0x0 C 



las cuales se pueden unir resultando, 

1 dE, DE, 



c dxo dxi c 



c 



(4.12) 



(4.13) 
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Las restantes ecuaciones de Maxwell son, 



dxi 



+ - 



C dXn 



O 

o, 



(4.14) 



quedando, 



dEk 1 dBi dBi 

OXj C OXo 



dxi 



d^J^^P = 0. 



(4.15) 



Finalmente las cuatro ecuaciones de Maxwell se pueden escribir en función 
del tensor electromagnético de la siguiente manera, 



c 



O 



(4.16) 



en donde se demuestra la ecuación vista anteriomente, DA" = — — con 
la ecuación de continuidad daA'^ = O y de ésta manera queda definida la 
covarianza del electromagnetismo. Con la covarianza de los campos electro- 
magnéticos, éstos transforman mediante las relaciones. 



E' = -f{E + l3x B) - 
B' = -f{B~ f3x E) 



7 + 1 

2 



7 



7 + 1 



/3(/3 . E) 
í3{f3.B). 



Problema: Sif<<cyf~c entonces qué pasa con 
las ecuaciones (14.171) ?: explique tanto matemática como 
físicamente. 



Problema: Si f = c en la ecuación (14.171) . qué significa 
para los campos electromagnéticos?. 



(4.17) 
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Capítulo 5 

Radiación Electromagnética 



El potencial escalar está definido por, 

^x, t)^ j j [t' + "ñZjñ. _ ^ (5.1) 

donde la contribución para el monopolo eléctrico se obtiene reemplazando 
|a; — x'l = = r, resultando, 



= ^ j d^x'p{x', t) J dt'S (t' + ^- í) 
= - j d^^'P t' = t--^ 



(5.2) 



pero al ser escalar, la contribución solo es esta. Para el potencial vectorial, 
bajo el gauge de Lorentz, está definido por (si desea desarrollar las ecuaciones 
en el sistema MKS o SI, deberá introducir un factor ^ ó según sea el 
caso). 



A 



(X, t)^\j d^x' J di' .S{t' + ^ - í) (5.3) 



donde el delta de Dirac hace que que los campos no se propagen instantánea- 
mente y por tanto se conviertan en campos retardados. Pero se observa que las 
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CAPÍTULO 5. RADIACIÓN ELECTROMAGNÉTICA 



fuentes varían en el tiempo sinusoidalmente de la forma p{x, t) = p{x)e 



-iüjt 



J{x, t) = J{x)e 
como, 



-iüjt 



y por tanto el potencial vectorial (Ec. (15.31) ) se modifica 



c 



d X dt 



.J{x')e 



\X — X' 



5[t' + 



\X — X 



- / d'^x 



3 , J{X') 



\X — X' 



dt'e-'^''6 i t' 



\X — X 



d^x'J{x') 



\X — X' 



d^x'Jix 



pik\x—x' I 



\X — X 



(5.4) 



Si la longitud de onda es A = y si las dimensiones de la fuente de radiación 
son mucho menores c¿ << A, entonces se definen tres zonas espaciales que 
definen la zona cercana o estática para d << r << A, la zona intermedia 
o de inducción para d << r ~ A, y la zona lejana o de radiación para 
d « X « r. Usando (ver Apéndice 19. 3p 



\x — X 



oo n 
4+1 



P;(C0S7) 



1=0 



°° ^ Ir'' 



(5.5) 



1=0 m=-l 



donde ^ = YlL=-i ^ *im {0\ (j)')YUe, 0) ([21], ¡23]), para la zona cercana 
kr << 1 y por tanto al analizar el límite cuando fcr — O, la exponencial 
tiende al valor unitario y así, usando (15. 5p y los datos del Apéndice 19. 3[ el 
potencial vectorial (15. 4p es. 



oo l 



iTT Y, 



Im \ 



^ 21 + 1 r 

1=0 m=-l 



l+l 



d^x'Jix'V'^Y * 



Im 



r, 0') 



(5.6) 



mostrando que en la zona cercana los campos son casi-estacionarios. 



Problema: Use el apéndice 19. 31 y demuestre la ecuación 

(¡ESD 



En la zona de radiación donde kr » 1, la exponencial de la ecuación (15. 5p 
oscila con gran magnitud, dando las características propias del potencial y de 
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los campos de radiación. Por tanto en el límite cuando kr ^ oo, se verifica 



r — ñ . x' = r 



(l — = r), Y por tanto el potencial vectorial (15. 4p es, 

1 r pik{r—ñ.x') 

A = - d'x'Jix" 



c J r — ñ .x' 

/ d3x'J(a;')e-*"-" 
c r J 



,r f:^?^/^Vj(.')(ñ.a.r (5.7) 

n=0 ■ 

lo anterior debido a que d « \ y por tanto se expande en términos de k. Si 
se tiene en cuenta solamente el primer término, el potencial vectorial es A = 
^ / d^x'J{x'), donde usando las propiedades de la derivada convectiva y de 
las soluciones fasoriales (^ = iu), que reescriben la ecuación de continuidad 
como V . J = iujp Y por tanto el potencal vectorial queda como A{x) = 
~^^^7~ / d^x'^'p{x') . Usando la definición del momento dipolar eléctrico p = 
qd = J d^x'x' p{x') y la definición del número de onda k = ^, el potencial 
finalmente toma la forma, 

pikr 

Aix) = -ikp . (5.8) 

r 

Finalmente los campos del dipolo eléctrico son, 

-ikp 

r 

= ]^{ñ^p)—[\-±-\ (5.9) 
r \ ikr ) 

y el campo eléctrico 

%c 

E = —V X B = ikV X (V X A) 

iü 

= (ñ X p) X ñ + [3ñ(ñ -p) — p] 



' e*'^^' (5.10) 
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ikr 



donde en la zona estática cuando se toma el lím^^^o y íisí ^y- — > O, (l — 
^, A;^ — > O y e*'^'' — ^) — ^(1 — ikr) ^, quedando los campos como, 



B = k'^{ñ X p) 



r kr 



ik , ^ \ 
= -(n X p) 



= — [3n{ñ . p) - p] , 



(5.11) 



mientras que en la zona de radiación cuando se toma el límfe^-^oo y Por tanto 
se conserva, pues decaen más o menos igual, mientras que ¿: ~^ O y 
~ ^) ~^ toman la forma. 



B 
E 



^ikr 

k'^iñ X p) 

r 

fc^(ñ X p) 



^ B X ñ 



X n 



(5.12) 



El vector de Poynting, que define la potencia radiada, está dado por S — 
^E X B*, pero lo que se necesita es calcular la potencia radiada per ángulo 
sólido unitario debida a la oscilación del momento dipolar se proyecta el 
vector de Poynting al vector unitario en dirección de r, ñ por que marca 
la característica del cambio de espacio, resultando. 



dP 
dñ 



2 

re 
k^c 



ñ.ExB* 



ñ .[{ñ X p) X ñ] X {ñ X p) 



Stt 
Stt 



[ñ X p) X ñ] 



h X psin^l 



pl^sin^e 



(5.13) 



donde resolviendo para el ángulo sólido se obtiene la potencia radiada total 



^ .S 1^1 • 



Capítulo 6 

Monopolos Magnéticos 



Si se analizan las ecuaciones de Maxwell (EMx) (ver ec. fl2.ll) ). se pue- 
de observar que éstas se dividen en macroscópicas (EMM) y microscópicas 
(EMm). Las EMm describen la electrodinámica que es válida para variaciones 
espacio-temporales arbitrarias del campo electromagnético acoplado, donde 
las fuentes describen todas las posibilidades; esto es: cargas puntuales, polari- 
zaciones, corrientes a nivel atómico, etc., en otras palabras describen el marco 
para distribuciones de fuentes y campos arbitrarios a escalas macroscópicas 
y microscópicas. Las EMM están definidas exclusivamente en cuerpos ma- 
croscópicos y por tanto para una correcta descripción electrodinámica de 
éstos, se definen unos nuevos campos llamados la densidad de fiujo eléctrico, 
o desplazamiento eléctrico D j la. intensidad de campo magnético o campo 
magnetizante H, siendo ambos originados de relaciones no-lineales, por lo 
general dificultuosas de analizar, de los campos E y B, las variables espa- 
cial, r, y el tiempo, t. En general se ha enseñado en los textos promedio que 
las relaciones entre las EMM y las EMm son tan sencillas como D = eE y 
H = fi^^B, pues el tratamiento de dichas cantidades se realiza en el vacío, 
lo que está muy por fuera de la realidad, pues en general éstas supuestas 
constantes llamadas permitividad e y permeabilidad /i se analizan dentro 
de materiales para visualizar las susceptibilidades eléctrica y magnética de 
éstos, resultando variables tensoriales de la forma. 



Se puede ver una falta de simetría en las EMx, que se puede solucionar 




(6.1) 
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si se suprimen las fuentes y por tanto la teoría electromagnética descrita se 
hace en ausencia de fuentes mediante las ecuaciones, 



V .E = 

V .B = 

V X E 



IdB 



c dt 
c ot 



O 



(6.2) 



donde la simetría buscada se logra pero con la consequencia de perder las 
fuentes, lo que es un problema grave pues si se quiere describir correctamen- 
te la naturaleza, las fuentes deben existir. Por tanto se debe establecer la 
simetría en presencia de fuentes, lo cual se intentó realizar casi desde la épo- 
ca de Maxwell, sin creer mucho en que las ecuaciones resultantes describirían 
algo en la naturaleza, sino que solo se verían simétricas matemáticamente. 
Posteriormente debido a Paul Dirac con sus famosos artículos del polo 

magnético |T], [10] y Cap. [TJ en los que definió al monopolo magnético como 
condicionante necesario para cuantizar a la carga eléctrica, con el inconve- 
niente de que resultaron divergencias diametralmente opuestas que definían 
una cuerda de divergencia, a la que se le dió el nombre de cuerda de Dirac, 
con la que algunos autores no están de acuerdo pues no se ha medido algo 
parecido en la naturaleza. Más recientemente, entre muchas investigaciones, 
algunos de los artículos que evidencian una dualidad electromagnética, ^12] . 
[T3] . [H], [15], [Tü], [TT], la relacionan desde el punto de vista topológico, aún 
teniendo en cuenta la cuerda de Dirac, e incluyendo teorías de gran unifica- 
ción como teorías de cuerdas o supercuerdas, y comportamientos no lineales 
como solitones o instantones. 

Para definir clásicamente al monopolo magnético, se realizan unas trans- 
formaciones entre los campos electromagnéticos (sin fuentes), cambiando B 
por —E y E por B, de la siguiente forma. 



V .E 

V .B 



O 
O 



VxE+-— 

c ot 

IdE 



V X B 



c dt 



E 
B 



B 



(V.B 
V .E 



E 



O 
O 



c dt 
IdB 

V X E + -— = O 

c ot 



(6.3) 
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Lo que no es natural, como se mencionó anteriormente. Si tal simetría no 
existiera entonces por ejemplo si hubieran habido dos campos Ei y Bi y se 
utilizaran otros campos E2 = B y B2 = —E, como una nueva teoría, al 
analizar la energía de las dos teorías, se tendría que 



Stt' ' Stt' ' 



1 

Stt' 



7^ ^1-^21^ + —|-B2|^ 

Stt 



y por tanto se verificaría que las cargas magnéticas no existen. Para mante- 
ner la dualidad electromagnética se deben verificar las transformaciones en 
presencia de fuentes, en donde se puede observar que E — > cB, cB — > —E, 



CPe 



ph y Ph 



-cpe, lo que marca transformaciones 



generales entre campos y fuentes electromagnéticas. Finalmente se tiene que 
dichas transformaciones llevan a las siguientes ecuaciones. 



V .E = Anpe 

V .B = ATiph 
^ ^ ^dB 

c ot 



4:71 

Jh 

c 



V xB- 



IdE _A'K 
c dt c ' 



E^B 
B^-E 

(pe, Je) 
{ph, Jh) 



'V.B = Anph 
V .E = -Aupe 



{ph, Jh) 

-{pe-, Je) 



Att 



— Je 

C 

An 

Jh 



^ ^dE 

c ot 

^ ^ ^dB 

c ot c 
(6.4) 

Se infiere que en el espacio 'normal' en el cual, en las ordenadas positivas 
está el campo E, en las ordenadas negativas el campo —E, en las abscisas 
positivas el campo S y en las abscisas negativas el campo —B, puede ser 
rotado una cantidad | en contra de las manecillas del reloj, dando como 
resultado el espacio dual en el que en las ordenadas positivas está el campo 
E, en las ordenadas negativas el campo —E, en las abscisas positivas el 
campo B y en las abscisas negativas el campo —B. Resumiendo éstos pasos, 
se verifican las transformaciones. 



B 



-E 



E = B 



(6.5) 



Igualmente se realiza está rotación a las fuentes electromagnéticas resultando 
las transformaciones, 

Pm = -Pe Pe = Ph- (6.6) 

Tradicionalmente se define el campo magnético mediante potencial vec- 
torial usando la relación. 



V . S = V . (V X A) = O, 
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lo cual es desafortunado, pues dicha ecuación niega al monopolo magnético, 
pues las líneas de campo ni convergen, ni divergen. Pero el potencial vectorial 
A es necesario para la descripción mecánico-cuántica de la electrodinámica, 
por tanto no se puede prescindir de él. Por tanto en un volumen r se debe 
verificar usando el teorema de la divergencia. 



V .BdT = A-ng = j B .dS (6.7) 

y por ésta razón el campo magnético debe estar definido por la ecuación 
S = V X A + X, donde x es un campo que evita que no existan los polos 
magnéticos y de ésta manera 



V . V X Adr + / V . = ^t^Q 
Jv 

V .xdr = A-ng. (6.8) 

Esta función, debe ser infinita en un punto y cero en cualquier otro lugar, 
para sostener la idea de una carga puntual en el espacio. Como podemos ver 
ésta es una característica típica de la función de distribución delta de Dirac 
(ver apéndice 19. ip . 

Si el volumen se escoge arbitrariamente, entonces x debe ser infinita en 
un punto de la superficie cerrada por el volumen formando una línea infinita 
que conecta al monopolo. Esta se llama la cuerda de Dirac, y debido a sus 
características, el campo S = V x A + % tendería a una singularidad, por lo 
que el potencial vectorial A debe compensar la singularidad en los alrededores 
de la cuerda de Dirac y debido a ésto no puede ser definido en cualquier punto 
del espacio. 

Hasta aquí todo parece correcto, pero el teorema de decomposición de 
Helmholtz dice que cualquier campo vectorial que sea lo suficientemente sua- 
ve y rápidamente decayente, se puede descomponer en un campo vectorial 
irrotacional (sin circulación) y en un campo vectorial solenoidal (sin diver- 
gencia) de la forma, 

T = -VT{V .T) + V xT(y xT) (6.9) 

donde se observa que si el campo es solenoidal, se recupera T = V x T(V x 
T) = V X Í7, y si el campo es irrotacional se obtiene T = — VT(V . T) = 
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Específicamente para la electrodinámica el campo magnético se escri- 
birá como, B = — V$ + V X A, lo que es válido solamente si V . -B y V x S 
tienden a cero más rápido que cuando r — oo, y si -B tiende a cero 
cuando r — > oo. De ésta manera A y $ son, 

1 f VxB . , 



J \x — x\ 

^ = — í ,^'^, c¿V (6.10) 

A-K \x- x'\ 



Problema: Explique por qué las ecuaciones ( ]6.10p son 
las formas correctas?. Es simplemente una solución ma- 
temática?. Existen más posibilidades físicamente admi- 
sibles? 



Entonces se definen V . S = Atx5^'^\x) y V x S = O pues si el campo 
magnético está definido mediante un polo, éste debe tener la forma ma- 
temática idéntica a la del polo eléctrico, B = g/ir~^f, y por tanto éste es 
irrotacional, lo que nos lleva a que el campo magnético esté definido por 
B = — V$ y por tanto para el polo magnético A = O, pues si existiera, vio- 
laría el teorema de descomposición de Helmholtz, ver ec. 16.91 y así aparece 
una línea de discontinuidad. Sinembargo el uso de 'V . B = Anph indica que 
no es correcto utilizar el potencial A. Adicionalmente la carga magnética se 
debe considerar como una distribución y por tanto la ecuación correcta debe 
ser V . S = ATTqh6^^\x). Si se escoge la existencia del potencial vectorial, 
entonces habrá una singularidad que no es necesariamente recta. 
Usando el teorema de Stokes, 



V .dL= a.dS (6.11) 

c Js 

donde a = V x V. Si se considera un flujo magnético, debido a un polo 
de la forma usual, entonces se observa que 4'{r,9) = JgB. dS = 27rg/¡(l — 
eos 6*) = §^A'.dL. Para que la definición del campo magnético, V x A', sea 
correcta sobre una superficie, ésta no puede ser cualquiera debido a la línea 
de singularidad. Así las posibles definiciones correctas del campo deben ser, 

^, Qh{^ — cosO) 7 

A = — ^ — -0 para 6 < n — e 

r sin 9 

A" = -^-h^±^4> Paraí^>. (6.12) 
r sin O 
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donde A' y A" están definidos de forma tal que evitan la singularidad de 
Dirac. A' alrededor del eje negativo de las z, y A" en cualquier punto del 
espacio, con excepción de 9 = e alrededor del eje positivo de las z. Los 
anteriores campos no serían potenciales vectoriales, sino campos vectoriales 
sin una definición física clara. El origen de estas ecuaciones se fundamenta 
en el análisis del campo de un solenoide largo y delgado. Ver Apéndice (19.51) . 
Con el objetivo de aclarar un poco más la teoría del monopolo magnético, se 
construye un campo magnético Coulombiano de la forma, 

B = g,^ (6.13) 

y por tanto una partícula con carga eléctrica qe en dicho campo sufre una 
fuerza magnética, Fm = QeV x B, de la forma, 

m-—r = — ^ -77 ^ (6.14) 



r3 \dt 

Suponiendo que la energía se conserva y solo es de movimiento, se puede 
proyectar la ec. fl6.14p en la trayectoria para analizar la órbita seguida por la 
carga eléctrica dentro del campo del polo magnético 

d'^r QpQh f dr \ q^qy. Id, , , , 

dt'^ mr^ \dt J mr^ 2 dt 

pero 

d'^r di dr\ í dr dr\ 1 d"^ ( '^'"'S^ 

= - „^ = o (6.16) 

donde se usaron los resultados anteriores para la última igualdad. Debido a 
que la potencia es nula se verifica la ecuación ^ (imf ^) = O que al resolverla 
(definiendo = 2Ct con C como la primera constante de integración y la 
segunda cero para í = 0) resulta la ecuación r = yvHF+W , que informa 
acerca del radio seguido por la partícula cargada, el cual no es cerrado en 
el sistema monopolar, sino que la cargas caen hacia el monopolo hasta una 
distancia mínima 6, y luego es refiejado hasta el infinito. 
Desde el punto de vista del momento angular L' = r x mv = mr x se 
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puede observar que 



dL d ( dr 

—— = m— r X — 
dt dt V dt 



r X m 



1^ 



r X 



QeQh 



fdr 
mr X — X r 
\dt 



QeQh 



QeQh ( dr 

dr\ 
r X m— X r 
dt ) 



qeQh 



L' X r) 

(6.17) 



donde se verifica que ^ . L' = O y por tanto ^ = O y que la magnitud del 



momento angular es L' = mvb. 
Si se analiza con un problema ordinario de Coulomb se observa que en éste 
problema el momento angular no es constante, como si lo es en un problema 
de tipo Coulomb. Modificando L' x r como, 



mr X 



dr 

Itt 



X r = m (r . r] 



mr 



dr _ dr 
3 I dt dt' 



dr 

d /r 



dr\ 
m \ r . — r 

dt ) 



dr dr 
mr I r— — r 

dt dt 



mr" 



dt \r 



(6.18) 



en la ec. f l6.17p . se llega a 

dL' q^qii ^ d fr^ dL 



-mr — — 



dt mr^ dt \rJ dt 

d / , r\ dL 

[L - qeQh 



dt 



r ) dt 



O 



y por tanto. 



cuya magnitud es 



L = L' - QeQh- = L' - qeqhf 



= L" - qlqi = {mvbY + {q^q.Y 



(6.19) 



(6.20) 



(6.21) 



resultando que la aparición de un término adicional en la definición del mo- 
mento angular genera una contribución no-trivial de campo. 
Como el vector de Poynting, S = -^E x B, tiene dimensiones de potencia 
por unidad de superficie, entonces el momento angular del campo eléctrico 
de una carga eléctrica puntual cuya posición definida por un radio vector r 
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y el campo magnético de un monopolo, está definido por, 



QeQh 



1 

Air 

qh_ 

Air 



d\-' [r' x{Ex B)] 



qh_ 

4:71 



r' X ¡Ex 



r./3 



E f' .E 



(6.22) 



Al integrar el primer término se observa que para un volumen que incluya 
el comportamiento de todo el sistema, ésta se anula. La segunda integral se 
integra por partes resultando. 



Finalmente usando V . E = A7rqe6'^{r — r'), se obtiene, 

Llq, = -qeqhr\ 



(6.23) 



(6.24) 



resultando que el valor del momento angular no es nulo aún en un sistema 
estático: carga eléctrica-monopolo magnético. 

Si se definen los campos como E = qei^ + E{qe2) y B = qni^ + B{qh2) 
para un sistema de pares de cargas (eléctrica- magnética) , llamado diones (d), 
entonces el momento angular se dicho sistema es. 



= ¿ / d'r' [r' X [E x B)] 



1 
1 

Air 



j dh' |r' X + E{q,2)) [qni^ + ^(5,2))] } 

J d^r' |r' X qei^ x B{qh2) + qhiE{qe2) ^ } 



(6.25) 



donde integrando de igual manera que en el caso anterior resulta un valor no 
nulo, L¿¿ = {qeiqh2 - qhiqe2) r. 

Si el sistema de cargas se acerca, éste sufre una dispersión del tipo Rutherford 
de la forma [22], [6], [9], 



da _ 1 
dñ^Q^ 



'2qeqh 



mv 



(6.26) 
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Explique por qué es una dispersión de ésta clase?. 
Compárela con la dispersión de Rutherford exacta. A 
cuáles conclusiones llega?. 



Todos los resultados anteriores muestran que el comportamiento físico de los 
monopolos es muy similar a lo que se observa en la naturaleza, y por tanto 
no es una incoherencia su descripción. 

Para analizar el potencial vectorial de un campo monopolar se usan con- 
sideraciones de simetría. Se supone que el campo magnético se comporta 
de igual manera a la ley de Coulomb, y por tanto la Lagrangiana del sis- 
tema está definida por L = ^mv'^ + q^v . A. Aplicándole las ecuaciones de 
Euler-Lagrange ^ ~ = O, se obtiene la ecuación de movimiento 



QeQh . 

— T^r X r, 



(6.27) 



donde el campo magnético se ha tomado como si fuera del tipo Coulombiano; 
ésto es, B = = 'V . A. Pero si se considera el campo magnético esférica- 
mente simétrico, el potencial vectorial correspondiente sería A = A{6)'V(j), 
donde usando la definición del campo vectorial de la ec. (16.121) . se elige la 
función A{9) como —qh{l + cos9), y por tanto después de calcular el gradiente 
Vy) = {r sin 6)~^ {— simpx + cos(fy) se obtiene el campo vectorial, 



qhjl + cosí 
r sin 6 



[siiaipx — eos Lpy); 



(6.28) 



en éste cálculo se ha usado la transformación de coordenadas esféricas, (p = 
— sin (fix+cos ífy. Finalmente se puede escribir el campo vectorial de la forma 



A{r) 



Qh r xn 



r r — r .n 



(6.29) 



donde el vector unitario n está dirigido hacia el eje z. Calculando el rotacional 
de éste campo se obtiene. 



B = V X A = V X 



Qh r xn 



r . n 



(6.30) 



lo que comprueba que el campo A obedece la forma de los campos de 
Coulomb o sea que es inversamente proporcional al cuadrado de la distancia 
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(r). También se observa que el carácter vectorial del campo vectorial desa- 
parece. Sinembargo para ^ = O o ^ = vr la ec. fl6.28p se anula o marca una 
singularidad. A ésto se le llama una línea semi-infinita de singularidad. Una 
posible solución a éste problema es reescribir el potencial vectorial como 



A(r) = (1 + cose)— U'^VU (6.31) 

Qe 

donde U = e~'^^'''^^'^ . De ésta forma el potencial es una transformación de 
gauge pura, la cual es singular y está complementada por un factor con 
dependencia polar. 

Analizando la vecindad de r, bajo un intervalo e {E? = + 6^), el poten- 
cial vectorial se puede escribir como, 

que en el límite cuando — > O, se obtiene el campo magnético 

B{r) = - AqhTiñe{z)5{x)5{y) (6.33) 
que resuelve todos los problemas con las ecuaciones de Maxwell (ec. (12.11) ). 

6.0.3. Teoría del Monopolo sin Cuerdas de Dirac 

Se esperaría que la teoría de monopolos y campos debería tomar una 
forma que es completamente dual a la teoría de cargas y campos y adicional- 
mente los campos electromagnéticos de un sistema de monopolos y los de un 
sistema de cargas deberían tener propiedades dinámicas idénticas, pero éstas 
propiedades deseables no se sostienen simultáneamente |18j, [19j . 

Las ecuaciones de Maxwell son = —AttJ^^^ y -^r- = O, con la fuerza 
de Lorentz dada por ma^^^-^ = q^Fj^^^v^, donde m es la masa en reposo de la 
partícula y Vy es su cuadri-velocidad. Todos los campos son derivados desde 
un cuadri-potencial regular mediante la ecuación, F(^e)^lu = ~ 
donde el cuadri-potencial A(^f.)a como generador de los campos cumple una 
importancia fundamental en la dinámica del sistema tal que forma la La- 
grangiana de interacción, C\-at = —J^^-^A(^e)^i- Intentando construir una teoría 



dual de la forma, = ~'^'^J[rn) y g^l"' = O, con la fuerza de Lorentz 
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"^^ím~) ~ l'^^í'm)'^'^- campos se derivan de cuadri-potenciales mediante 



"(m) ~ yní-' (m) 
{iri)^iv fe'' dx" 

forma, (ver Cap. [5]), A(e) = —iAe^'^^~~^^i, los campos están definidos por 



^tnúav ~ a "/^^" Tfir^- Teniendo en cuenta un potencial vectorial de la 



E=-^ = cüAe'<^-'n (6.34) 
dt ^ ^ 





B = V X A(e) 



= ujAe''^y--'>j. (6.35) 

Con el monopolo magnético el potencial dual se define de igual forma, como 

^(m) = — '¿v4e*'^^~~*)j, y aplicando las transformaciones duales E B y 
B —>■ —E, se hallan los campos. 




= cjAe'"(^-*)j. (6.36) 



B = -?^=uAe'<^-')i. (6.37) 
Por tanto se oberva que 

E = ^ = -V X A(^) = E (6.38) 

y que 

B = Vx = -^B, (6.39) 

concluyendo que se pueden identificar los campos de radiación de cargas con 
los campos de radiación de monopolos, lo que hace innecesaria la definición 
de los monopolos con cuerdas de Dirac. Se puede observar adicionalmente 
que debido a que la acción es un escalar de Lorentz, todos los términos de 
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la densidad Lagrangiana lo deben ser también; como la electrodinámica es 
lineal, el término de interacción carga-monopolo debe ser una suma de can- 
tidades bilineales que contienen dos factores: uno relacionado con cargas y el 
otro con monopolos y que un sistema de una carga sin movimiento - mono- 
polo sin movimiento, no cambia en el tiempo. 

Las cargas no interactúan con campos ligados de monopolos y los monopo- 
los no interactúan con campos ligados de cargas. Las cargas interactúan con 
todos los campos de radiación de las cargas y con los campos de radiación de 
los monopolos. Los monopolos interactúan con todos los campos de radiación 
de los monopolos y con los campos de radiación de las cargas. Los campos de 
radiación de las cargas y de los monopolos son considerados como una sola 
entidad y se denotan con el subíndice w para formar el tensor ligado de cam- 
pos de radiación de cargas, Fl^^^-^ y el tensor ligado de campos de radiación 
de monopolos, F^^^y Finalmente la forma de la fuerza de Lorentz ejercida 
sobre las cargas es, ma^^^ = 'imF¡^^^^V(^e)v Y la ejercida correspondientemente 
sobre los monopolos es, rna^^^^^ = <imF^^^-^V{m)u, donde a los campos mono- 
polares se les denomina campos magnetoeléctricos [3], [1]. 
Para desarrollar una teoría de monopolos de ésta clase, se debe tener en 
cuenta que ésta debe reducir a la teoría conocida de cargas y ondas para un 
sistema sin monopolos y a una teoría dual para sistemas donde las cargas 
no existen. La teoría se debe originar de una densidad Lagrangiana cuyos 
términos son funciones regulares de los potenciales, de sus derivadas, de las 
corrientes de carga y de monopolos. 

Los términos de dicha teoría deben ser escalares de Lorentz y deberían serían 
invariantes bajo traslaciones espacio-temporales asociadas al grupo de Poin- 
caré; adicionalmente las ecuaciones de movimiento de sus campos deben ser 
lineales y deben sostener el principio de superposición. 
Por último se debe asumir que un sistema de una carga y un monopolo, 
donde las dos partículas no se mueven, no cambia con el tiempo. 



Capítulo 7 



La Teoría de los Polos 
Magnéticos 

En éste capítulo se pretende explicar y desarrollar detalladamente las 
ecuaciones y los procesos matemáticos del artículo original de P. A. M. Dirac 
[10]. Lo que se busca aquí es presentar la teoría original del monopolo y su 
cuantización con el objetivo que el lector entienda los conceptos originales y 
la matemática involucrada que llevaron a Dirac a la regla de cuantización de 
la carga eléctrica, usando la teoría de polos magnéticos. 
El artículo de Dirac [TÜ]: 'La teoría de los polos magnéticos', el autor describe 
la forma general de las ecuaciones clásicas de movimiento enmarcada en el 
tensor electromagnético y construye la dinámica electromagnética mediante 
la ecuación de Lorentz, así como la teoría de los polos magnéticos compa- 
rando continuamente con el caso eléctrico de forma tal que conceptualmente 
adiciona términos a las ecuaciones para solucionar problemas relacionados 
principalmente con las cuerdas de discontinuidad o cuerdas de Dirac, sin que 
éstos afecten las ecuaciones de movimiento. El método usado es el de cambiar 
el punto de referencia al de las líneas de mundo de las partículas y nuevas va- 
riables sin carácter físico que orientan dichas cuerdas contenidas en láminas 
asociadas a cada polo. Los potenciales retardados se definen en éstas láminas 
pues son contribuciones de los polos y son del tipo Lienard-Wiechert. El autor 
modifica tales potenciales para que sean consistentes con la teoría básica. El 
paso siguiente es la definición del principio de acción de la electrodinámica, 
que debe ser modificado para que no hallar inconsistencias con las cuerdas 
de Dirac usando la condición que una cuerda nunca debe pasar a través 
de una partícula cargada. Esta modificación da como resultado un término 
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relacionado con la variación del potencial que no afecta las ecuaciones de 
movimiento. Con todo lo anterior el autor define la función Hamiltoniana 
para poder cuantizar la teoría y por tanto se escriben los momentos tanto de 
las partículas cargadas como de los polos usando la técnica de Fourier. Las 
ecuaciones de movimiento están descritas entre dos superficies tridimensio- 
nales: de los campos y de los polos, de forma tal que las partículas no existen 
sino hasta cuando llegan a la superfice de los campos y lo mismo para los 
polos suponiendo que éstos están unidos a las cuerdas y que su superficie 
tridimensional vas má alia. Sinembargo a pesar que algunos campos dejan 
de evolucionar dinámicamente, el resto del sistema continua hasta que todo 
éste evoluciona totalmente. Este comportamiento para las acciones de ciertas 
partículas o campos y otros no es útil para el manejo de la teoría pues, el 
sistema no tiene que detenerse dinámicamente cuando solamente un campo 
lo hace. Con las anteriores consideraciones se construye las ecuaciones que 
relacionan la dinámica de las coordenadas y los momentos (tanto para las 
partículas cargadas como para los polos) y con éstas se hallan las ecuaciones 
de Hamilton-Jacobi con las que se cuantiza la teoría, obteniendo la condición 
de cuantización de Dirac que define la necesidad de la existencia del polo 
magnético para poder cuantizar la carga eléctrica. 

I-II. Introducción y ecuaciones clásicas de movimiento 

El autor describe la elecrodinámica desde el punto de vista de la simetría 
de los campos eléctricos y magnéticos concluyendo que para que ésta sea 
coherente el polo magnético debe exisitir y el soporte de dicha afirmación es 
que hay teorías exitosas que se han basado en partículas predichas teórica- 
mente y que después de cierto tiempo se han descubierto. Un ejemplo para 
las que se descubrieron son los neutrinos y para las que no, es el bosón de 
Higgs. También se apoya en el hecho de que el polo magnético (monopolo) 
es necesario para la cuantización de la carga eléctrica, dada por, 

eg = ^nhc (7.1) 

donde n es un entero y (7 es la carga magnética; por tanto la carga eléctrica 
para su cuantización requiere que sea definida mediante 
En la sección Las ecuaciones clásicas de movimiento se hace una breve des- 
cripción de la teoría de los campos electromagnéticos usando tanto el ten- 
sor electromagnético F^j, = — F¡^^ (ver Capítulo Hj), como su dual Fjj, (En 
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el Capítulo m el tensor dual está definido por JF'^^), con las propiedades 
Fo\ = F23, FJ3 = — Foi y F^l = —F^y. Adicionalmente se usa un nuevo 
tensor electromagnético G¡j_y que cumple la condición F^^G^'^ = F^yG^^^ . 
Las ecuaciones de Maxwell están dadas por = —Anj^, donde es la 
cuadri-densidad de corriente eléctrica y para el espacio dual se cumple que 

la cuadri-divergencia se anula; esto es = O, resultando una falta de si- 
metría en las ecuaciones de Maxwell que es corregida si la ecuación en el 

espacio dual se define como = — 47rfc^, donde kf^ es la cuadri-densidad 
de corriente magnética. Usando la línea de mundo de la partícula cargada 
que determina su evolución espacio-temporal y está definida por = Zfj_{s), 
donde s es el tiempo propio de la partícula, una carga puntual e que se des- 
plaza en una línea de mundo infinita genera una cuadri-densidad de corriente 
eléctrica dada por 

e 

Para el caso magnético, la cuadri-corriente magnética está definida de forma 
similar: 

k,{x) = Y.3 j ^^'i^ - ^)ds. (7.3) 

Las funciones delta 5^{x) = S{xo)S{xi)6{x2)6{x3) (ver apéndice 19. ip son usa- 
das para describir un punto de discontinuidad en el espacio de las cargas y 
describe la suma sobre las partículas cargadas (e) o sobre los polos (g). 
Con el fin de describir la dinámica del sistema se escribe la fuerza de Lorentz 
F = e{E + V X B) en función del tensor electromagnético quedando. 

Dicha ecuación describe el comportamiento dinámico de una partícula eléctri- 
ca cargada a lo largo de su línea de mundo y la descripción de los campos 
se realiza en el punto z donde está situada la partícula. Para una partícula 
cargada magnéticamente, usando consideraciones de simetría, se asume una 
cuadri-fuerza de Lorentz, 



(Pz^ dz" 



m^ = g—FUz). (7.5) 
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Debido al carácter de los campos si la medición de éstos se hace exactamente 
en el punto z, dicha magnitud originará una divergencia; por éste motivo las 
variaciones se deben realizar en magnitudes bajas, pero debido a que la fun- 
ción delta varía brúscamente, es conveniente reemplazarla por una función 
más suave que cumpla las propiedades de la función delta o imponer con- 
diciones límite tales que la teoría solo es invariante de Lorentz en el límite 
mismo y después de éste, pero no antes (ésta es la elección hecha por Di- 
rac). Dirac modificó el tensor electromagnético F^,^ al tensor F*^ que cumple 
las condiciones requeridas, resultando las ecuaciones de movimiento para las 
cargas eléctricas y para los polos. 



47r^e y ^6\x-z)ds (7.6) 

e 

4^E^/ ^S\x-z)ds. (7.7) 



dF* 
dxy 

dF'f ^ r dz 



Las cuadri-fuerzas de Lorentz también se modifican, para las partículas car- 
gadas y para los polos a, 

d'^z dz"" 
d'^7 dy^ 

donde solo falta estructurar al nuevo tensor electromagnético. 



III. Potenciales electromagnéticos 



Para crear una teoría que pueda llevarse a la mecánica cuántica, se nece- 
sita llegar al principio de acción (ver apéndice 19.61) mediante las ecuaciones 
de movimiento y para ésto se necesitan los potenciales electromagnéticos. 
Mediante la definición F^¡^y = — en función del cuadri-potencial, se 

dF dF^ y 

observa que se cumplen las ecuaciones = —^i^j^ y -q^ = O, pero éstas 
no cumplen con la condición del polo magnético pues su cuadri-divergencia es 
nula. Por dicho motivo la definición del tensor electromagnético debe variar 
de alguna manera debido a que éste falla en algunos puntos de la superfi- 
cie que encierra al polo, formando una 'especie' de cuerda que se extendería 
desde el polo hacia el inifinito uniendo todos los puntos sobre cada una de 
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las superficies que encierran al polo; el inconveniente es que dicha cuerda 
no tendría significado físico y debido a que es infinita no existirían variables 
para describirla. 

Las variables necesarias para fijar las posiciones de las cuerdas son coordena- 
das dinámicas y momentos conjugados a ellas. Dichas variables son necesarias 
para una teoría dinámica pero no corresponden con observables y no afectan 
el fenómeno físico. 

Cada cuerda esquemaliza una lámina bidimensional en el espacio-tiempo; 
éstas láminas serán las regiones donde el tensor electomagnético F^¡j falla y 
cada lámina podría ser descrita expresando un punto general y^, sobre ella 
como función de dos parámetros tq y ri (y^ = ?/^(ro, ri)). Se supone que cada 
lámina se extiende al infinito y por tanto los parámetros tq y ri pueden ser 
asignados con valores de forma tal que ri = O sobre la línea de mundo del 
polo y se extiende al infinito siguiendo la cuerda y — oo < Tq < +00 como 
valores que describen desde el pasado infinito hasta el futuro infinito. 

La característica de ésta cuerda es que en sus extremos deben estar loca- 
lizados los polos. Así el nuevo tensor se define como, 



una de las láminas y su suma se realiza sobre todas las láminas cada una 
de las cuales está asociada a un polo. Como se mencionó anteriormente, 
cada una de estas láminas que encierran al polo puede ser descrita por un 
conjunto de coordenadas (ro,ri), formando una función = y^{To,Ti) con 
el objetivo de orientar los puntos en los cuales existe la divergencia debido a 
la cuerda. Reemplazando la ecuación fl7.10p en la ecuación (17.71) . se obtiene 
el siguiente desarrollo: primero calculando la cuadri-divergencia del tensor 
electromagnético del polo 




(7.10) 



9 



donde el campo Gj^^ es un campo que desaparece en todo sitio excepto sobre 




9 




(7.11) 



9 



y por tanto. 




(7.12) 



a 
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pero debido a que Gj^^ = —G^i, y a que d'^Fj^^ — —4:nJ2g9 J ~ 

entonces, —^T^Yl/g^^^iJ-f — ~^'^Ylig9 1 ^'^^(^ ~ ^)^^ J P^^ tanto es obvio 
que 

d'G,. = g J ^S\x - z)ds; (7.13) 
su solución se logra usando el teorema de Stokes 

\ox oy oy ox J J \ ox oy J 



para x ^ tq, y ^ r^, U ^ S^{x -y)yV^y^,, 
^d6^{x-y) dy^ dy^, dS"^ {x - y) 




drQ dri dri Otq 



(ÍTqTi 



reorganizando el primer miembro se llega a 

' dy^ d5^ {x - y) dS^ {x - y) dy^ 




dri dro dro dri 

s 

usando la regla de la cadena para y¡y 

' dy^ dyy dy^ dy^ \ d6^{x -y) 



+ ^^dn] (7.15) 



(ÍTqTi C^-IG) 




dn dro Oto dn ) dy^ 



dTon (7.17) 



5 

ahora se intercambia la variable de integración para la función delta quedando 

y por tanto queda 

d f f fdy, dy. _ dy^ dy,\ _ ^^^^^^^ 



dx¡yj J \(9ri dro dro dri 
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Finalmente al incluir el segundo miembro del teorema de Stokes y teniendo 
en cuenta la carga magnética se obtiene, 



dx 



V 



9 j S\x - y) (^dro + ^dn^ (7.20) 



recordando que í/^ = y^{To,Ti) y tomando ri = O, debido a que se integra 
solo sobre una sola lámina del anillo y que no tienda al infinito, entonces 



I S\x - y) (^4^) dro = d^G,^, (7.21) 



haciendo las correspondencias tq = s y = ^/¿(to^O) es la línea 

de mundo de la partícula, se observa que corresponde con la ecuación ( 17.13^ . 
También aquí se puede observar que el campo que está definido solamente 
sobre las láminas y que le da consistencia a la electrodinámica conel polo 
magnético incluido tiene la forma: 

Las ecuaciones (17. 6p y (17.71) estn directamente relacionadas mediante la ecua- 
ción Ffj_iy = 1^ — + Att J2g Gj^j, para definir los potenciales retardados 
(PR) en las líneas de mundo de las partículas previamente definidas y su con- 
tribución a cada partícula depende solamente de la línea de mundo de ella y 
de la lámina que está unida a ésta en el caso que halla polo. La contribución 
de los PR se define mediante una función invariante de Lorentz definida por 

J{x) = {^^^^"^ paraxo>0 ^^^3) 
O para Xq < O 



o por 



J{x) = -6{xo-r), (7.24) 



r 

r = [xl + xl + xl)^. (7.25) 

La función A{x) de Jordán y Pauli (ver apéndice 19.20 se relaciona con J{x) 
mediante 

A{x) = J{x) - J{~x). (7.26) 



52 



CAPÍTULO 7. LA TEORÍA DE LOS POLOS MAGNÉTICOS 



Se verifica directamente alrededor del origen usando la integral de nj{x) 
sobre un pequeño volumen cuadridimensional alrededor del origen como una 
integral superficial tridimensional sobre la frontera de ese volumen que, 

nj(x) ^Anó^x). (7.27) 

La contribución de una partícula cargada a los PR está definida por, 

K,rM-^eJ J{x-z)^ds (7.28) 

basada en los campos y potenciales dados por Liénard-Wiechert, 

A''{x) = ^ J d^x'D,^t{x - x') J"(x'), (7.29) 

donde Dj.f¡i{x — x') es el propagador retardado (ret) o avanzado (av), 

l^^^et ^e{xo-x'^)5[ix~x'f] ^^^^^^ 
[¿ = av ^9{x'q - xo)S[{x - x'f], 

y donde 9{xo — x'q) cs la función de Hcaviside o función escalón. La contri- 
bución correspondiente de un polo está dada por la integración sobre toda la 
lámina, y debido a que su origen tiene un carácter antisimétrico, debe estar 
presente el tensor totalmente antisimétrico Sap'ys- 

A í\ [ [ dy^ dyp dJ{x -y) 

A.U^) = ge.,,. J J ^^^^^^oár, (7.31) 

Para evidenciar como es la variación de ese PR en la lámina, se realiza una 
variación cuadridimensional quedando, 

s'-^A.M^) = // f f ^W^'^^"'*^" 

Usando la relación 
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la ecuación anterior se modifica a, 

dy^ dy"^ d'^J{x — y) dy^ dy^ d"^ J [x — y) 



(9ro dTi dyi^dy/s Dtq dn dyi^dya 
dyP dy'^ d'^J{x-y) dy^dy^d'^J{x-y) 



dro dri dyt^dyp dro dri dyi^dy^ 

dya Qyt^ Q2j(^^_ Qya QyP Q2j [x-y) \ ^^^^^^ 



dro dri dyi^dyp ' dro dn dy^'dyp 



rr f gy" d'^Jjx - y) 9/ d'^Jjx - y) dy^dy^^j^ 
^JJ [dri drody^ dr^ dndy^ dro dn ^ ^ 
dy^d^Jix-y) dy-d^Jix-y) dy<^ ^y^.^j^^ _ ^ ^^^^^^ ^^ ^^^ 



dri drodya dro dndy^ dro dn 

f [ ¡dy^d'^J{x-y) dyf^d^J{x-y) dy^d'^J{x-y) 
^ jj \dTi drodyg dro drodya dri drodya 
dy^^Jix-y)\ ffídy^dy" 



Qya QyP 



/3 ^ 

-Ó4{x -y)> drodri, (7.36) 



dro dri 

donde □J(a; — y) = ^ti5^{x — y). Por tanto usando, JJ ^ ^dlodyf drodri — 
I dy"-'-^ y dy-'-^ - í^l '-£ds, se obtiene 



J y=z 



_dJ(pv) , _ + íí {^^S.(^ - y) 

dya y 13 J J J y dro dn 



Qy» QyP ^ 

~^"^<^4(^ - y) I drodri. (7.37) 
Integrando mediante el uso de la función delta y debido a que, 

J l dya y/3 J 
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es un intercambio entre a y ¡3 y teniendo en cuenta que, 

J { dro Oti dro dn 

la anterior ecuación se escribe abreviadamente de la forma 



dy^ dy°' dy"' dy^ i 

5i{x — y) — — 5i{x — y)\ dr^dri (7.39) 



9 



^-^^^dy- + ^^^r'^dy^ 



dyp 



dya 



al3 



(7.40) 

cambiando la referencia a la línea de mundo de la partícula descrita con las 
cuadri-coordenadas Zfj,{s), donde s es el tiempo propio de ésta y e^'^°'^dfiAi,^^ct{x) 



la ecuación se transforma a 



ret 



dJ{x — y) 
dJ{x — z) dz' 



dz"" 



dz. 



ds 



■ds — (a/?). 



ds — (aP) 



(7.41) 



Usando el hecho que = 9 J 



dJ{x-z) dz°' 



dz, 



■0 



ds 



ds, O sea que se mantiene la corres- 



pondencia del potencial de Lienard-Wiechert con los potenciales magnéticos 
para un campo retardado producido por un polo magnético (campos magnéti- 
cos retardados -B^t); ecuación finalmente toma la forma 



ret 



dXB 



+ 



dB 

dXr 



ret 



(7.42) 



IV. Principio de acción 

El principio de acción (ver apéndice 19.61) se enuncia sobre tres términos: 
el de la partícula (Ji), el del campo {I2) y el de las interacciones de las car- 
gas eléctrica y magnética con el campo (Js). Entonces la acción se define 
I = I1 + I2 + h, donde 



Ji = m J ds 



e+g 
1 



/2 = — y F;.(a:)F^^(x)d^x 
h = Y.^¡ A'i^)^ds (7.43) 
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donde = J Fi^i,{x')^{x — x')d'^x' con 7(x — x') definida como una función 
par que se aproxima a la función delta 5^{x) y que tiene como objetivo 
evitar infinitos en la ecuación de movimiento debido a los campos infinitos 
producidos por las cargas puntuales y por los monopolos magnéticos; por ésta 
razón el delta I2 se redefine como I'i — SS Fniy{x')^{x — x')F'^^ {x')d'^xd'^x' . 
La variación de la acción 7i está dada por, 



5h 



d'^Z 



ds^ 



^Sz^Ss. 



(7.44) 



Para I3, 



5/3 = ^^ 



dx^^ dx'^ 



+ SA^ 



dz'^ 
ds 



ds 



(7.45) 



Finalmente la variación en la acción del campo es 



5/2 = [F^A^)6F^''{x)^{x - x') + 5F^,{x)F^^{x)-í{x - x')] d^xd^x' 

= ¿ / [f;^{x)5F'^''{x)] d'x (7.46) 



pero Fi'^{x) = - |^ + 47r Gt/^'^ y 5F'''' = 
por tanto reemplazando de obtiene 



dxu dxv 



+ 47rE,<^Gt'^^ 



1 

47r 



5h = -— / F;^{x)-:—d'x 



dx^ 



X 



1 r dF* (x) 1 r 

4^ J ^^^^"'^''^ + 2 5^ y Fl:{x)SG,.d'x (7.47) 



Usando la definición Gfj,y{x) = 2g JJ ^^^4{x — y)dTodTi su variación 
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está dada por, 



SG^^ix) -'¿gjjs (^^) ¿4(x - y)dTodn 



, o 1 1 (^y^^ ^^4(2: -y) 



(7.48) 



por tanto SI2 es 



57^ 



1 

47r 



+ 




dSy.dy^ , dy^dSy 



Oto Oti 

dy^dyf,d64{x - y) 
Oto Bti dy" 



h{x - y) 



Sy" 



drodrAéx. (7.49) 



Integrando sobre x, 



51' = — 



1 

47r 



^x)— — d X 
ax,j 



Y.- 



FtM 



dro dri 

dy,dy,dFl:{y) 



d5y^dy^, ^ dy^,d5y^^ 



-Sy" 



dTodTi, (7.50) 
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pero i-MMSy^] = '-^Sy^ + FlM^, entonces, 



dro 



SIo 



— í F* ( 



,4 , V- 

X}— — a X + 9 

9 



dx^ 



dóy^dy^ , dyf^dóy^ 



+ 



dyadya ( d 



dro dri dro dri 
dóy^ 



d 



dSy^ 



o 



dro dn dro ^'^i'-^y'^y J aro ''"^ ' dr^ 



drodn (7.51) 



Siguiendo el desarrollo algebraico, se realiza el cambio a la línea de mundo 
tal como se hizo anteriormente y usando el teorema de Stokes, 

dU dV dU dV 



^ ddro ddri ddv 
se obtiene la acción total. 



' dV \ f fdV dV \ 



(7.53) 



donde. 



dxt^ dx'^ 



dz^ 
ds 



ds 



dFt (9Ft 



+ 



vp 



dyp dyi^ 



dF^^t , dF^; dFl; \ dyP dy" 



dyP dy^^ dy'-" J dro dri 



-óy'^drodT, = O, entonces ^ + ^ 4 



t* 



O, ésto significa que todos los puntos se mantendrán en la lámina. 
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Igualando a cero los otros coeficientes se puede observar que se reproducen 
las ecuaciones de movimiento iniciales, lo cual le da consistencia a la teoría 
y se concluye que las cuerdas de Dirac nunca deben pasar por las partículas 
cargadas. 

De = O se observa que al aplicar el principio de acción no se llega a una 
ecuación de movimiento lo que significa que las cuerdas de Dirac no poseen 
ecuación de movimiento y por tanto no tienen carácter físico. 
Con la acción anterior existe el problema de la anulación del momento con- 
jugado a Aq en la formulación Hamiltoniana. La solución de ésto se logra 
adicionando un término a la acción de la forma, I4 = J ^^^^d'^x, cuya 
variación es 



'4 

1 f dAt dSA 



í^d^x, (7.55) 



47r J dxn dxi, 

integrando por partes con U = dV = ^^^dx'^, dU = -^^¡^^dx'^ y 

V = óAf^ y usando la condición de integral extremal, se obtiene. 



pero debido a la condición suplementaria que |^ = O o |^ = O, entonces 

el término „ no afecta las ecuaciones de movimiento. En la formulación 
Hamiltoniana es necesario distiguir entre las ecuaciones que se mantienen 
en virtud de las condiciones suplementarias y las que son independientes de 
éstas; entonces. 



dxj, 



-47r^e j ^S4{x-z)ds (7.57) 



con F,. = 1^ - 1^ + in^^Gl y G,. = g JJ (^^ - |^^) S,{x - 
z)dTodTi entonces reemplazando. 



dxydx- 



^ + 47r y ^-4^ye [^54x- z)ds, (7.58) 

^ dx^ ^ J ds ^ ' ^ ' 
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debido anula por la condición suplementaria. Finalmente se 

obtiene el D'Alembertiano, 



0A; = A^J2'J ^6,ix - z)ds + A^Y: (7-59) 



V. El método para pascir a la formulación Hamiltoniana 

Para cuantizar una teoría se usa convertir las ecuaciones de movimiento de 
un sistema dinámico en forma de principio de acción en forma Halmiltoniana. 
El método es tomar la integral de acción previa a un tiempo t y hallar su 
variación permitiendo la variación de t. Entonces la variación de la acción 
SI es función lineal de 5t y Sq en las coordenadas dinámicas del tiempo. Los 
términos adicionales en SI se cancelan cuando se usan las ecuaciones de mo- 
vimiento. Se introduce la variación total en los q finales, Aq = 6q + qót y 
usando la transformación de Legendre 61 = J2Pr{Qr,Qr)A.qr — W{qr,qr)St, 
se definen los momentos Pr y la energía W que son función de las coorde- 
nadas generalizadas qr y de las velocidades generalizadas g^- Estas variables 
son constantes en el tiempo si la Hamiltoniana transformada es nula y que 
satisface la relación 

W - H{p,q) = (7.60) 

que es la ecuación de Hamilton-Jacobi cuyo argumento son derivadas parcia- 
les de la acción respecto a g y a í. Debido a que existen varias ecuaciones que 
relacionan a g y a p, entonces deben existir varias ecuaciones de Hamilton- 
Jacobi que originen varias ecuaciones de onda. 

Relativistamente es el mismo procedimiento pero se evalúa la integral de la 
acción sobre el espacio-tiempo relacionado con una superficie tridimensional 
que se extiende al infinito. Lo anterior se puede modificar debido a que en 
lugar de detener la integral de acción en un tiempo determinado o en una su- 
perficie S espacial-oide tridimensional definida, se podrían detener diferentes 
términos en la acción en diferentes tiempos. Lo anteriormente expuesto es 
posible si se describe la dinámica del sistema enmarcada en difentcs partes a 
diferentes tiempos y después de que algunas partes del sistema han dejado de 
evolucionar voluntaria o involuntariamente desde el punto de vista dinámico, 
el resto del sistema lo sigue haciendo hasta que eventualmente, termine la 
dinámica del sistema. Cada vez que una parte del sistema para, se puede de- 
finir para ésta una ecuación de Hamilton-Jacobi. Este método no siempre es 
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conveniente para solución de muchos problemas por tanto se generalizará pa- 
sando a un espacio en el que una superficie espacial-oide tridimensional S 
en la que se puede detener diferentes términos de la integral y no toda su 
acción. El hecho es que si se evalúa cuidadosamente la dinámica del sistema, 
se puede observar si parte de éste se comporta de manera no natural o no 
física; una vez hecha la evaluación se detiene esa parte de la integral de acción 
y puede hacerse dicho proceso a diferentes tiempos de evolución del sistema 
y el resto de éste seguirá su desaroUo regido por la ecuación de movimiento 
hasta cuando ésta se detenga, colapse, etc. 

Una forma de detener la integral de acción cuando existe interacción entre 
las partículas y los campos es suponer que las partículas no existen en el 
espacio-tiempo en donde se traslapan los conos de luz de diferentes partícu- 
las y posteriormente detener el campo donde el intervalo de tiempo para 
hacer ésto debe ser de gran magnitud para dejar que éste se estabilize. En- 
tonces se varía la integral de acción detenida mediante la variación sobre los 
puntos de la línea de mundo de la partícula, 2;^, en el espacio-tiempo donde 
ésta desaparece y también en la superficie tridimensional de los campos Sp, 
donde el campo deja de existir. Intentando emular el principio de Hamilton 
o el teorema fundamental del cálculo variacional, se iguala a cero la parte 
de la variación de la integral de acción detenida y que no está conectada 
con las variaciones de frontera, se obtiene la ecuación de movimiento de las 
partículas antes que desaparezcan y de ésta manera se originan las ecuaciones 
de campo que rigen a la partícula aún después de haber dejado de existir. 
Debido a las variaciones realizadas en las líneas de mundo, Az^^, localizadas 
dentro del campo, se obtienen ecuaciones más apropiadas para manipular 
que las del método tradicional en donde las partículas y el campo dejan de 
existir simultáneamente. 

Con la nueva electrodinámica se supone que todas las partículas y también 
las cuerdas unidas a los polos dejan de existir en la superficie tridimensional 
de los polos Sp, y que los campos electromagnéticos dejan de existir mucho 
más tarde sobre la superficie de los campos Sp, resultando que las integrales 
de acción dadas por las ecuaciones Ji e I3 de la ec. 07.431] . tienen que dete- 
nerse cuando las líneas de mundo alcancen a Sp, mientras que las acciones 
dadas por las ecuaciones I2 de la ec. (17.431) y (17.561) . dejan de existir en la 
frontera Sp. El hecho de que se detengan éstas cantidades no significa que el 
sistema colapse, pues el resto de la dinámica de éste seguirá evolucionando 
hasta Sp y aún después cuando llegue a Sp. 

Si la conexión entre U y U* dada por U*{x) = j U{x')^{x — x')d^x', y es tal 
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que el valor de cualquiera de ellos en un punto x es determinado por los valo- 
res del otro sobre los puntos espacio-temporales cercanos a x formando una 
especie de recursión. Así si uno se anula en cierta región espacio-temporal, 
el otro también lo hará en esa región excepto en las cercanías de ésta. 
Debido a que G^^ se anula en cualquier lugar excepto sobre las láminas, 
debe hacerlo también en la región entre Sp y Sp-, excepto en los puntos cerca 
de donde las cuerdas dejan de existir. En ésta región se tiene que deja de exis- 
tir la primera suma sobre DA* = 47r e J — z)ds + Att 
debido a que la integral se detiene en Sp; por tanto DA*^{x) =0, y se según 
lo analizado también se tiene que □y4^(a;) = O en la región entre Sp y Sp 
con excepción de los puntos cercanos a donde las partículas cargadas dejan 
de existir. En las regiones donde los D'AIcmbcrtianos de y A*^{x) se 
mantienen se pueden definir éstas mediante transformaciones de Fourier de 
la forma: 



á* 



feo 



donde kx = koXQ — kiXi — k2X2 — k^Xs, d^k = dkidk2dk^, k^ = ±(A;^-|-A;|-|-A;|)^ 
y la suma se realiza sobre los valores positivos y negativos de k^. El factor 
^ se introduce debido a que éste es invariante de Lorentz. La condición que 
y A* sean reales da las condiciones, 

A_,^ = -^^, A\^ = Al^. (7.63) 

Sea la resolución de Fourier de la función ^{x) dada por, 

l{x)^^j^ié'^dH (7-^ = 70, (7.64) 

entonces la condición de función par ^{—x) — ^{x) da 7_; = 7; debido a que 
7; e M. Se tiene entonces que 
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donde usando el delta de Dirac se obtiene, 

.d^k 



feo 



ko 



ikx _ 
Jkx'^ 



O 



^k.e''''^ = (7.66) 



Es necesario conservar la transformación de Fourier de ^^(x) en cada punto 
z donde una partícula cargada desaparece y mantener la transformación de 
Fourier de A*^{x) en y donde la cuerda deja de existir. Es probable que lo 
anterior se pueda manejar de tal forma que al elegir una función 7 que origine 
una condición para que el punto y nunca esté demasiado cerca a debido a 
que ésto significaría que la totalidad del sistema deja de existir. Si se asume 
un campo U{x), esté definido por U*{x') en los puntos x' muy cercanos a los 
puntos X y por fuera del cono de luz de x. Así la transformación de Fourier 
de A^{z) es definida por la transformación de Fourier de A*^{x') sobre los 
puntos x' para los cuales la transformación de Fourier de A*^{z) es válida 
y está definida. De la misma forma la transformación de Fourier de A*^{y) 
será válida solamente si U*{x) está definida por U{x') en los puntos x' muy 
cercanos a a x fuera del cono de luz de de x. 

La condición suplementaria = O, = O, se modifica en la región entre 
Sp Y Sp- Con las integrales dadas en las ecuaciones (I7.28P y (17.311) detenidas 
en Sp y usando z' para z{s'): 

dAl^ret _ r dJjx-z') dzl^^, 



dJ{x — z') dzl^ , 



I dzl ds' 

= - ^eJ{x-z) (7.67) 

e 

donde — eJ{x — z) es diferente de cero debido cuando x esté sobre el cono 
de luz futura de cualquier punto z donde una partícula cargada no exista. 

VI. La formulación Hamiltoniana 



Dejando variar a, Sp j no a, Sp con el objetivo de formar la variación de 
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la integral de acción acotada justamente como en la sección anterior y eva- 
luando los términos en 51 conectados con los límites, resulta que los términos 
originados en la electrodinámica, 



e+g ^ ' e 



(7.68) 



donde /S.z^^ son los cambios totales de las coordenadas del punto donde la 
partícula deja de existir y que vienen de ÍJi y (5J2. Ya no es posible usar la 
acción fl7.47p sino la acción 



„ J — oo 



X 



X 



(7.69) 

La razón del cambio en el segundo término es que si Sp está lo suficientemente 
lejos de Sp, entonces 7(0: — x') = O para los x anteriores a ¿"p y los x' 
posteriores a. Sp- Se usarán los métodos que se usaron para obtener fl7.54p 
sobre las láminas que se extienden solamente sobre las partes de las láminas 
anteriores a Sp, y entonces se tendrán extra-términos, en forma de integrales 
de línea a lo largo de las líneas donde se encuentran las láminas, originados 
debido a la aplicación del teorema de Stokes. Organizando la parametrización 
de las láminas de forma tal que la línea donde la lámina encuentre a Sp 
es dada para tq = constante y la línea donde la lámina variada halle la Sp 
variada es dada para tq = la misma constante, éstas integrales de línea toman 
la forma 



9 



dri 



(7.70) 



Las líneas de integración son las posiciones de las cuerdas cuando dejan de 
existir. Para la definición de ÍJ4, ya no es posible usar (17.551) . En su lugar se 
usa 



6L 



1 id Al 85 A^" dA^dóA^"* 
^TT J-00 V (^^^ 9xi, dx^^ dxy 



d X, 



(7.71) 
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que es parecido al primer término de 



SI' = — 



Stt 



F: 



^SA^' 



d X 



dx, ' dx, 

+ i E / {Pt^^G^"' + F.^^G^'"*) d'x; (7.72) 



sumándolos da 
Stt 



1 





- F* 1 







asAf" 

dx¡, 



+ 



dA, 
dxi^ 



^5A^'* 
dx„ 



éx + 



(7.73) 



Usando F^ 
cantidades anteriores como 



dAl 35 A^^ _ dAl dAl ^ dA*^ 
^x^' dx„ ^x^' dxi^ dx" 



dA* 



y también 

dA^_p d5Ai^_dA^_dA^^dA^ 
dxf^ dxi, dxi^ dx^^ dx^ 



dA,^ 



(7.74) 



+ - V Gj,, = - 47r V Gt 



Por tanto la ecuación queda 



1 

Stt 





[( 


1 —oo 





dx" 



^5A^' 
dxy 



dA 



dx 



^5A^'* 
dx^ 



(7.75) 



d^x + 



1 

Stt 





[( 


1 —oo 





dx- dx, l^^i-^ 



dx, j 



+ 



9 



teniendo en cuenta solamente los potenciales, entonces 



1 

Stt 



3o \ dx'' dx, J 



d^x + 



ídA^dSA''* 
3o V dx'^ dx. 



d X 



+ 



dS + - 
(7.76) 

(7.77) 
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integrando por partes el primer término usando U = -g^, dV — 



d X, 



entonces dU 



dx°'dx- 



-^^d^x, V = 5A^^ y por tanto, 



dx" dxy J 



d'^x 



¡9 A* 1 
^ "dAf" 
ox" 


Sp 
— oo 


\dAl 1 
^ ^ÓA^" 
ox" 


Sf 

— oo 


/ ^ ^ÓA^dS 
J dx" 



ÓAf"- — J^d'^x 

-oo 



dx°'dx 



V 



Sp g2^* 

" ÓAf'dS'^dx" 



dx'^dx 



V 



(7.78) 



donde se ha integrado sobre la superficie tridimensional 5*^ usando d^x = 
dS^dx" y se ha usado el hecho que los potenciales se anulan en el infinito y 
en Sp (condición de integral extremal). El resultado es 



1 



57r 



57r 



dA^ 

"ÓA^ + ^5A^* 1 dS"" 



dx'' 



dx'^ 



57r 



. 98 A^' 85 A^'* 



dx„ 



fJ,U 



dx„ 



dA* dA„ 

^'SA^ + —^ÓA^"* dS" 



dx'' 



dx'^ 



dóA^" . ^6A^'* 
h G' 



d^x + 



d^x + 



(7.79) 



Los restantes términos en 61, se cancelan debido cuando se usa la ecuación 
de movimiento, sin que Sp esté muy cerca a Sp. Por tanto se llega a. 



+ 



Stt 



g 

dA* dA 

^SAf' + ^SA^'* ]dSr (7.80) 



dx'^ 



dx'^ 



Pero con la expresión anterior no es posible encontrar los momentos con la 
transformación original 61 = 'Y^p.^^Qr — W6t, pues las variaciones de A'^ y 
A^* dadas en fl7.79p no son independientes entre si y debido a que la superficie 
de los campos Sp no ha variado. 
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Para hallar los momentos se pasa a los componentes de Fourier de los 
potenciales para los que se pudiera utilizar las transformaciones de Fourier 
dadas por las ecuaciones fl7.6ip y fl7.62p . debido a que existe interés en la 
expresión fl7.79p con los potenciales sobre la superficie Sp- Se considera un 
movimiento tal que satisface las ecuaciones de movimiento de forma tal que 
las transformaciones de Fourier (17.611) y (17.621) son válidas sobre Sp también 
para dicho movimiento variado. Así la expresión (I7.79P se convierte, hallando 
las variaciones de los potenciales, 



fco 

ik'x d?k 



pero debido a que 5A^^ = J2k¡j I ^^k'^^^^'^^^ entonces, 

^-^5A^ = iy íí KAl SA''^e'('+''> — ^ 

fco,fc¿ 



d^k d^k 



usando A^ = 'jkAki,, entonces 



^^SA- = fc„A»,7..í^"-e'<-«->í^í|: (7,83) 



fco,fcó 

y su suma es 

?^5A^ + ^5A^* = z y /"/" fc,(7fc + 7fc')^fc SA^'^e''^'^''^—^. 
dx'' dx-" ^ jj ^Kík^ ikj k, ^ ^, 

fco,fc¿ 

(7.84) 

Al llevar la integral a la superficie tridimensional para 5"^ tal que xq = 
constante, se tiene que el último término de la ecuación (17.801) se modifica a 



Sn^lll ^^^^^ + ^'^')^'^M''^^'^'^''^''^^dS'^ (7.85) 

fco,fc¿ ^ 
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y por tanto solucionado parcialmente ésta ecuación se llega a 

ítt^ ^ // (^-t + lk')Ak^5A'''>^e'^^'+^'^^''^5\k + k')d^k^. (7.86) 

Este resultado se logra con xq ~ constante y k,^^^ — hko^^ — kidS'^dxi, 
para Xi se integra resultando, 

/e*.-V^.,„,.,.2.,(,„ + ,i>...*i.» (7.87) 

El resto de las integrales para X2 y xs tienen resultados idénticos matemáti- 
camente, resultando finalmente 

pues 6%k) = 6{ki)6{k2)6{k3). 

Se observa claramente que el delta de Dirac hace que el integrando se 
anule para todos los puntos excepto para k¡ — —ki, lo cual implica que 
k'o — ±/co- Entonces usando 7_; = 7; la integral se reduce a 



T íí Í^^^T^) ^ik,M'Ve*('=o+fcí,)xo¿3(^ ^ k')d'kd'k' (7.89) 

con 7_i = 7¿, si /c¿ = —ki y /cg = ±ko entonces 7^/ = 7-^ — y por tanto 
¿^(/c + A;') = 1. Así se integra dando 

^^"E / (z|) ^^.'^^''"^'^ (7-90) 

y si si A;q = +A;^ y A;^ = — A;¿ entonces se integra dando 

= ítt^ I /7^^1W,\ _,,e^^'=°-d^^, (7.91) 



,2 



feo,fc¿ ' 
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finalmente la solución es 



= -2z^ y: i (^) A,Mk'.d'k 

+ (7-92) 

feo 

Pero (7fe + 7feo,ifci) ^feM<^^feo -fef = ^ (Tfe^fe^^feo -fef) > P^es 7^ + -fkoM = (^7fe J 
Ak^SAi^ _f,i:^ — SAk^, lo que forma un diferencial exacto quedando 

feo feo 

(7.93) 

El segundo término se descarta porque es un diferencial exacto y no afecta 
las ecuaciones de movimiento. Finalmente, 



J (z^) (^•94) 

A_k^ = -Ak, 

Si ko > O, usando A*_j^^^ = —A^.^ entonces la ecuación se modifica a, 

7-1 = li 



= -2z7r2 5^ |7fe {AuM'"- - ^fe/^'') (7-95) 
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Debido a que A^^^ {^i^^SA^^) ^-A+fc^7+fcM ^ 7-fc^-fc„ = IkAk^^A^^ y 
usando 8 {A^^A^^) = 6Ak^A''^' + AkJA^^^ se llega a 6 [Ak^A^^) -2Ak^5A^^ = 
ÓAk^A'^- -aI.ÓA'^'^ Y as( 

= -2^n'6 j IkA.^A'^^ - j ^^KM'^'^^ (7.96) 

donde el primer término de la integral no se tiene en cuenta debido a las 
razones anteriormente expuestas. Finalmente la variación de la acción queda, 

= E + E ^Mz)^z' ^^^r (^) "^^^ 

-AtTc' I^.A.^ó'^^, (7.97) 



donde son las coordenadas dinámicas de las partículas cuando dejan de 
existir; y^{Ti) son las coordenadas de los puntos sobre las cuerdas cuando 
éstas dejan de existir y que construyen un continuo unidimensional de coor- 
denadas para cada polo; para cada valor de /i y A^^ son las componentes 
de Fourier (con > 0) de los potenciales después que las partículas y las 
cuerdas han dejado de existir. Los coeficientes en la ecuación fl7.97p serán 
los momentos conjugados. Así los momentos de una partícula cargada son 



m- 



ds 



Pfj, = rn-^ + eAij_{z) y los de una partícula con polo serán 
El momento conjugado a las variables de la cuerda (3^{ti) = y^{ri) son 
í^^iji) = Q^lv^ y ^1 momento conjugado a está dado por —^^'^^Ak^. 
El momento í3^{ti) forma un continuo unidimensional de variables que corres- 
ponde al continuo unidimensional de las coordenadas y^{Ti) y el momento de 
campo conjugado a forma un continuo tridimensional que corresponde 
al continuo tridimensional de las coordenadas de campo. 

Para las coordenadas y los momentos de cada partícula se definen los 
brackets de Poisson de la forma 

{Pm, Zu} = l'rn^ + eA^,{z), zA = \m^, zA + {eA^{z), z^} 



ds \ \ ds 

dz^ 
ds 



rn^,Zy\ = g^,y, (7.98) 
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para las coordenadas y los momentos de la cuerda 

{P''in),y,ÍT[)} = g^Jin-r[), (7.99) 
para las variables de campo 

|.— . [ iin^ — ^ iko f . 

{Ak„Aki\ = < — —lkAk^^,Aki \ = {A^^^, Ak'J 

iko 3 



g,,5\k - k') (7.100) 



y todos los otros brackets de Poisson son cero. 

Usando — eA^(z) = Y Pfi = ""^^5 con el fin de eliminar la 

velocidad para cada partícula cargada, se obtiene 

\p, - eA.izf = [p, - eA.iz)] [p^ - eA^{z)] = (m^^ = 

[p^-eA^{z)][p^ -eA^'{z)]-m^ = (7.101) 
y para cada polo la ecuación es 

p^p^' = m^ ^ p^p^ - = O, (7.102) 
que deben estar unidas con 

í^'{n)~9Fl:{y)^^=0, (7.103) 

donde los potenciales Afj^{z) en las anteriores ecuaciones están dados por las 
componentes de Fourier, Ak^{z) y Ak^{z). Las ecuaciones (17.1011) . (17.1021) y 
(I7.103P solo relacionan coordenadas y momentos dinámicos. Son ecuaciones 
diferenciales satisfechas por la integral de acción I cuando los momentos se 
ven como derivadas de / y son ecuaciones de Hamilton-Jacobi de la teoría 
de los polos magnéticos. Las condiciones suplementarias (I7.67P deben ser 
tratadas como ecuaciones de Hamilton-Jacobi y éstas ecuaciones obtenidas 
tomando diferentes puntos del campo x no son independientes de las ecuacio- 
nes de movimiento o entre si y de esta manera se obtiene una serie completa 
e independiente de ellas haciendo una transformación de Fourier en la región 
entre Sp y Sp- En ésta región, usando la función de Pauli- Jordán (I7.26p . se 
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pede reemplazar J{x — z) por A{x — z) cuya trasformada de Fourier está dada 
por 

ko 

entonces se obtienen, 

e e 

que son las ecuaciones que envuelven solamente variables dinámicas de coor- 
denadas y momentos y que son conformadas de la manera correcta para 
construir las ecuaciones de Hamilton-Jacobi. 

VII. Cuantización 

Para cuantizar se reemplazan las coordenadas dinámicas y el momento de 
la teoría clásica por operadores que satisfagan las reglas de conmutación 

{P'{n),y,ÍT[)} = g,J{n-T[) 

{Ak,:A,,} = - k') (7.106) 

y entonces se reemplaza las ecuaciones de Hamilton-Jacobi por ecuaciones 
de onda obtenidas igualando a cero las miembros izquierdos de las ecuacio- 
nes de Hamilton-Jacobi, que están escritas en función de operadores , y que 
que están aplicadas a la función de onda í^Ij. Las ecuaciones de onda obte- 
nidas de ésta manera son consistentes entre si debido a que los operadores 
sobre ip en sus miembros izquierdos conmutan y esto es debido a la anula- 
ción de los miembros izquierdos de los brackets de Poisson de las ecuaciones 
de Hamilton-Jacobi. La solución da como resultado ecuaciones de onda para 
partículas sin cspín; para los electrones se reemplazan por las ecuaciones de 
onda de cspín Debido a la falta de información sobre los polos, se supone 
con espín ^h. 

Con las matrices de espín usuales ai, «2, «3, «m; para cada partícula car- 
gada existe una ecuación de la forma, 

{po - eAo{z) - ar\pr - eAr{z)] - amm} ip = O, (7.107) 
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para cada partícula con un polo, 

[po - «rPr - amm] ip = O, (7.108) 

para cada cuerda, 

ip = 0, (7.109) 
y para las variables de campo, 

(^^7.^,. - ¿ 5^ ee-^'^^ j ^ = O, (^^7.1,. - ¿, ee^^'^^ ^ = 0. 

(7.110) 

La función de onda i¡j se podría tomar como como una función de las variables 
de la partícula z^, de las variables de espín disponibles para cada partícula, 
de las variables de la cuerda y^{Ti) con rango O < ri < oo y de las variables 
de campo A^^. Dicha función estará definida solamente cuando los puntos 
^A" VÁ'^^) 6stén fuera de los conos de luz de las otras partículas. 
La ecuación (17.1081) sugiere que el campo electromagnético no actúa sobre 
los polos pero si sobre las cuerdas (ec. (17.1091) ) y debido a que los polos están 
restringidos a los extremos de las cuerdas entonces el campo no afecta el 
movimiento de los polos. 

VIII. La carga y el polo unitarios 

La integral de acción de la teoría clásica podría ser considerada como una 
función de puntos espacio-temporales donde las partículas dejan de existir, 
de las líneas y^iri) en el espacio-tiempo donde las cuerdas dejan de existir 
y de las variables de campo apropiadas. Esta acción está definida siempre 
que las cuerdas no pasen a través de cualquiera de los puntos z^ donde las 
partículas cargadas dejan de existir. 

Se mantendrán todos los puntos de las partículas {z^) fijos y también fi- 
jas todas las cuerdas excepto una que variará continuamente manteniéndola 
siempre en la superficie tridimensional Sp situada alrededor de uno de los 
puntos z^ donde una partícula cargada está situada justamente antes que 
deje de existir y luego regrese a su posición original. Al mismo tiempo los 
potenciales A^{x) son variados continuamente con el objetivo de preservar 
las ecuaciones (I7.10p y (I7.22p con valores fijos para el campo F^y(x) y que re- 
toman los valores originales junto con la cuerda. Debido a que la deformación 



/3.(n)-í7F¿(y)|- 
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no puede ser restituida continuamente en la acción entonces una cuerda no 
puede pasar a través de una partícula cargada. La cuerda se extiende sobre 
una superficie cerrada bidimensional a, permaneciendo en Sp y encerrando 
el punto donde se sitúa la carga y dicha superficie no puede ser encogida 
continuamente a cero debido a que no debe pasar a través de una carga. 
Sinembargo se espera variar la acción I bajo éste proceso de deformación. 
Una pequeña variación de una cuerda y de los potenciales con los fijos 
lleva a una variación de I {DI) , dada por la suma de los miembros derechos 
de las ecuaciones fl7.7ip y fl7.72p dada por 



DI = — 

Stt 



Sf 



dx,. 



+ 



dxf^ 



dx,. 



d*x 



„ J — oo 



(7.111) 



debido a que los tensores electromagnéticos normal y dual se mantienen fi- 
jos y que los cuadri-potenciales normal y dual son regresados a sus valores 
originales, entonces la variación de la acción se modifica a 

DI = - Y^ {F^^JG^^" + F^JG^^''*) d'x, (7.112) 

que como ya anteriormente se explicó se puede reparametrizar como (ver la 
ec. (E70])) 

'dy" 



DI 



E 



9 



FllSy' 



dri 



ÚTi. 



(7.113) 



Se observa que para un proceso de deformación cerrado los parámetros de 
los que depende el punto sobre la lámina ya no son importantes debido a 
que no se desea la descripción de la línea de mundo del polo y por tanto 



DI = g ¡ FllSy^dy^ 



(7.114) 



donde óy'^dy'^ es un elemento diferencial de la superficie bidimensional exten- 
dida por la cuerda y estará definido por da'^'^. De esta manera la variación 
de la acción será, 



DI 



9 



(7.115) 
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La integral resultante tiene la forma matemática de una divergencia y por 
tanto definirá el fiujo eléctrico total que pasa a través de la superficie a dando 
como resultado: DI = Auge. 

Se podría rodear cualquier carga con una cuerda por un número indefinido 
de ciclos y por tanto la incertidumbre total en / está dada por la suma (suma 
sobre todas las cargas e y sobre todos los polos g) J2ge ''^ge9^, donde rrigi, 
es un coeficiente arbitrario que define un análogo al peso estadístico de cada 
uno de los términos de la suma. 

Asociando el comportamiento no univaluado de la integral de acción / de 
la teoría de los polos con el de la acción de algunas de la cantidades físicas, 
se puede observar que es posible usar la regla de quantización de Bohr para 
cada carga y polo, resultando (para c = 1) 

4nge = nh, (7.116) 

donde n es un entero. Usando el momento angular de Planck h = ^ \a 
ecuación queda como 

n 

ge = -hc. (7.117) 

Este resultado puede ser obtenido por la teoría desarrollada en la capítulo 
VII. Cuantización sin la regla de cuantización de Bohr que comprende 
cuatro postulados que determinan: la forma de la órbita del electrón alrededor 
del núcleo (sin radiar energía electromagnética) y su relación con la Ley 
de Coulomb, su limitación entre la infinita cantidad de órbitas (mecánica 
clásica) y solo en las que el momento angular es un múltiplo entero de la 
constante de Planck y el cambio de órbita mediante la radiación de energía 
electromagnética. 

Aprovechando el segundo postulado se puede usar la condición que la función 
de onda debe ser univaluada. Para las coordenadas y momentos de la cuerda, 
cuyo conmutador está dado por la ec. (17.991) . {/5'^(ri), í/í/(t{)} = g^u^i^Ti —t[), 
(usando el hecho que I3^{ti) = ih g^ ^^^^^ , debido a que éste es un operador 
momento) y donde se observa que el momento de la cuerda puede ser definido 



mediante la ec. dTTTÜÜ]) , P^in) - gFll{y)£-^ 



= O como 



mostrando la manera como la función de onda varía cuando las coordenadas 
de la cuerda se modifican. Si una cuerda es desplazada y se curva sobre una 
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superficie bidimensional cr, usando la ec. ( 17.118P se puede observar que 

f = l^^'-'^'f ■■■ / f = -f /^-'^'Ir*^ 

Integrando se obtiene 

por tanto la función de onda ip queda multiplicada por 

g-f /Ftíd..^^ (7.121) 

donde se puede observar que F^* ocurriendo en diferentes puntos del integran- 
do conmutan entre si, lo que podría ser una analogía a una transformación 
local de campos. Debido a que ip es univaluada, éste debe retomar su valor 
inicial y por tanto la fase (17.1211) debe ser unitaria. Para que esto sea posible 
se debe cumplir la igualdad | / F^^da^" = 27m siendo n un entero. Se puede 
observar que ésta condición vuelve a definir la ecuación (17.1161) . Ange = nh. 
Finalmente la conclusión a la que se llega es que [lOj: 

...la cuantización de las ecuaciones de movimiento de las partícu- 
las cargadas y de las partículas con polo, es posible solamente 
condicionando a las cargas y a los polos a que deben ser múlti- 
plos integrales de una unidad de carga cq y de una unidad de polo 
qq satisfaciendo, e^go = ^hc... 



(7.120) 



CAPÍTULO 7. LA TEORÍA DE LOS POLOS MAGNÉTICOS 



Capítulo 8 



Incoherencias en la Teoría del 
Monopolo Magnético 



Durante muchos años los físicos han intentado hallar el monopolo in- 
fructuosamente. El intento de explicar las discontinuidades inherentes en la 
teoría de Dirac o cuerdas de Dirac y su resultados negativos en la medición 
experimental, han llevado a artificios físicos y matemáticos que han hecho 
que la teoría sea poco natural. Actualmente existen sistemas de monopolos y 
campos que son duales a la electrodinámica clásica (EDC) con los cuales se 
pueden explicar la supuesta falta de simetría de ésta (ver capítulo (16.0.31) ). 
A pesar de ésto muchos científicos han seguido investigando y usando las 
cuerdas de Dirac en la teoría del monopolo. 

Las ecuaciones de Maxwell, que gobiernan el movimiento de los campos 
electromagnéticos son. 



y en ellas están contenidas las ecuaciones de Maxwell y la relatividad especial 
de Einstein (ver Capítulos [3] y H]). La fuerza ejercida sobre toda la materia 
que posee carga está definida por la fuerza de Lorentz en cuatro dimensiones 
o la cuadri-fuerza de Lorentz, 



donde el tensor electromagnético F^'^, está definido según lo dicho en el 



d^F^ = - 
d^F*/" = O 



(8.1) 



QeFrv, = ma>¿ 



(8.2) 
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Capítulo m El tensor electromagnético dual (ver Cap. HI) se define como 

F*^"' = ^e'"''"^F^p (8.3) 

el cual es totalmente antisimétrico, justamente como F^'^ también lo es. La 
relación intrínseca entre estos dos tensores y la definición de la densidad 
Lagrangiana 

^ -FrF,,^, - J^A,,^ (8.4) 



I6n 

lleva a realizar transformaciones duales que conducen a las ecuaciones de 
Maxwell duales, las cuales son usadas para definir la teoría de un sistema de 
monopolos con sus campos electromagnéticos de la forma, 

-a.Fr = O, (8.5) 

en las cuales están contenidas las ecuaciones de movimiento de Maxwell para 
los campos monopolares. La cuadri-fuerza sobre la materia monopolar es de 
forma análoga a la ejercida sobre la materia cargada 

qHF:t''Vu = ma^^. (8.6) 
La densidad Lagrangiana del monopolo está dada por 



IGtt 



Con éstas dos teorías se puede observar que [E] ninguna de éstas relaciona 
las dos cargas (eléctrica y magnética), pues cada teoría contiene las ecuacio- 
nes de movimiento así como la ecuación de la dinámica para una carga y no 
la relaciona con la otra. Dichas teorías deberían conformar una nueva teoría 
carga eléctrica-monopolo magnético en la cual, en los límites en los que una 
de las cargas no exista sobreviva la otra. Pero en la teoría de los monopolos 
de Dirac si se calcula el límite en el que la carga eléctrica no exista, no se 
recuperan las ecuaciones de Maxwell en el espacio normal. A pesar que la 
teoría de Dirac se construya bajo la suposición (deseable) que tanto la carga 
eléctrica como el monopolo tengan propiedades dinámicas iguales, ésto lleva 
a inconsistencias físicas tales como las cuerdas de Dirac. 
Por ejemplo los campos electromagnéticos radian con las condiciones E'^ = 
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B'^ Y E . B = O que son bien conocidas en las ondas electromagnéticas y 
cuyo significado es que los campos se mueven oscilatoriamente de manera 
perpendicularmente a la velocidad de radiación o de propagación de la onda 
electromagnética lo que muestra de nuevo un grado de simetría que está en 
total desacuerdo con una singularidad de cuerda. 

Adicionalmente aunque con la teoría de Dirac al modificar a 'mano' algu- 
nas de las ecuaciones de la dinámica electromagnética y/o de la acción (ver 
el capítulo [7¡) para reproducir la electrodinámica llamándola nueva electro- 
dinámica y que no es posible derivarla de una función Lagrangiana elec- 
tromagnética corriente por lo que las variables canónicas no estarán bién 
definidas y por tanto la función Hamiltoniana será muy difícil de obtener; 
adicionalmente la carga magnética transforma como un seudoescalar y no 
como un escalar lo que dificulta el entendimiento físico. 

Existen muchos artículos de investigación donde se describen los incon- 
venientes de los polos de Dirac para casos específicos como el de Daniel 
Zwazinger [2Ü] donde explica que éstos son prohibidos en la teoría de la ma- 
triz de dispersión debido a que dicha teoría es invariante de Lorentz y los 
polos de Dirac no respetan dicha invarianza lo cual está muy bién determi- 
nado en [5] donde analiza la no-covarianza del monopolo de Dirac resultando 
una incompatibilidad con la teoría cuántica de campos que para solucionarla 
habría que redefinir muchos conceptos básicos de ésta. En fin hay mucha lite- 
ratura al respecto pero el problema fundamental es la definición física de las 
cuerdas de Dirac y el entendimiento matemático de la aplicación del teorema 
de Stokes debido a que éste no se define en todo el espacio por la divergencia 
generada por la cuerda de Dirac. 
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Capítulo 9 
Apéndices 



9.1. Delta de Dirac 

La función delta de Dirac es una función de distribución que permite de- 
finir puntos matemáticos de discontinuidades en el espacio. Un ejemplo muy 
claro es su uso para definir una carga puntual. El delta de Dirac está de- 
finido por las relaciones 6{x) = O para x = O y / 6{x)dx = 1 si x 7^ O y 
donde la región de integración incluye x = 0. Para una función arbitraria 
/ f{x)6{x)dx = /(O), donde la región de integración incluye x = 0. La repre- 
sentación del delta de Dirac está dado por 5{y) = lím/i^oo ES]' [¡25]. 
La integración se logra usando 

r e(— = lím e(— = 2 lím ^^^¡^(^^ " ^-)] 

= 2nS{{cü,-a,)). (9.1) 

La relación de ortonormalidad se construye 

j u*{r)uj{r)(fr = 6{u^ - a^)6{ujy - cry)6{uj^ - a^) = 6^{u - a). (9.2) 

En un intervalo {—x,x) siendo a una cantidad positiva se impone la condi- 
ción, 

(1 a a 
— para < x < — 
\ 2 (9.3) 
O para |x| > — 
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de forma tal que en dicho intervalo la función delta tomará el valor ^ for- 
ma una especie de pico. Si se integra sobre las x en ese mismo intervalo 
y sobre una función arbitraria f{x) (bien definida en x = 0), entonces se 
tiene que J_^^ S°'{x)f{x)dx y teniendo en cuenta que f{x) ^ /(O) entonces 

J_^¡^ 6°'{x)f{x)dx ^ /(O) f^^ S°'{x)dx = /(O). Si a es muy pequeña entonces 
la aproximación será mejor y por tanto al evaluar el límite cuando tiende a 
cero se verifica 



6{x)f{x)dx = f{0), (9.4) 

J+x 

y de forma más general 

P — OD 

5{x)f{x-Xo)dx = f{xo), (9.5) 



' +00 

Algunas de las propiedades de la función delta de Dirac son: 

6{x) = 6{—x) 
d6{x) d6{—x) 

djCC 

x6{x) = O 
d6{x) 

X ; = —ÓÍX) 

dx 

6(ax) = -5{x) 
a 

5(x^ - a^) = ^[^(a; - a) + 6{x + a)] 



6{a — x)6{x — b)dx = 6{a — h) 

f{x)6{x — a) = f{a)6{x — a) (9.6) 
Usando la ecuación (19. 5p y 

F{k) = (27r)-5 / e-'^''F{x)dx 



oo 

+ 00 



F{x) = (27r)-3 / e+"'''F{k)dk, (9.7) 

J — oo 

la relación entre el delta de Dirac y Fourier está dada por 

/ + 00 
e-*^^5(x - Xo)dx. (9.8) 
-oo 
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En los textos l2il. 1261. 125 



se puede profundizar más sobre éste tema. 



9.2. Función de Pauli- Jordán 

La función de conmutación de Pauli- Jordán (FPJ) está dada por, 

Aix-y) = t[ipix),ip{y)] (9.9) 

donde 



A(x) 



(2vr)= 



dk 



r (9.10) 



donde se verifica que su variación temporal está relacionada con la función 
delta de Dirac de la forma: 



aA(x°,x) 



5(x). 



(9.11) 



9.3. Ecuación de Laplace en Coordenadas Esféri- 
cas (CE) 



La ecuación de Laplace en CE se escribe como, 



1 



d ( . ^dT\ 1 d'^T ^ 
smé»— 1 + „ . „ „ ^r-Tr = O 



^ ¿^^2 ^2 gjj^ 06 \ 06 / sin^ 6 ■^'^^ 



(9.12) 



Si se asume una solución de variables separables de la forma J-' = ^^^0(^)$(0), 
al reemplazarla en la ecuación, y después de reordenar los términos ésta que- 
da, 



sin^ 6 



1 d'^R 



+ 



1 



d f . je 

I si]2 6' 

R dr^ ' r2 sin ^9 d6 \ d6 



1 



0. 



(9.13) 



Si se elige una constante de separación — m^, entonces se verifica la ecuación 
diferencial, 

(9.14) 
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cuya solución es, 

* = e^'""^ (9.15) 

donde m debe ser un entero, pues si no lo es, la función $(0) puede ser 
multivaluada. El resto de la ecuación es entonces, 



r^sin^^ 



1 d^R 



1 



d 



R dr"^ sin 6*6 d9 



~ —m 



y por tanto 



d'R 



+ 



1 



R dr"^ sin ^0 dO 



smi 



^dQ 



m 



(9.16) 



(9.17) 



Separando variables de nuevo con una nueva constante de separación /(/ + 1), 
se llega a 

r'^d^R 

+ (9.18) 



despejando se obtiene. 



R dr^ 

d'^R l{l + 1) 



R 



(9.19) 



donde l{l + 1) es otra constante. La solución de ésta ecuación está dada por 
funciones hiperbólicas, 

i? = 7^le^+^+7^2e-'. 



finalmente se tiene la última ecuación diferencial. 



d 



sinOede 
la cual se reorganiza a 



m 



1 d 
sin 6 d6 



sin 9 



dQ 

Ib 



-l{l + 1) + - 



m 



e. 



Si se realiza el cambio de variable x = eos 6, entonces dx = 
— (1 — x^) 2(i^, y sin^ 9 — 1 — x'^. Por tanto la ecuación queda, 

i2 



d_ 

dx 



'l-x' 



d& 

dx 



-1(1 + 1) + 



m 



x^ 



9. 



Finalmente haciendo P{x) = 0(x), se tiene la ecuación. 



d_ 
dx 



-/(/ 



m 
1-x^ 



(9.20) 
(9.21) 

(9.22) 
sin 9d9 ~ 

(9.23) 
(9.24) 
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a la que se le llama ecuación generalizada de Legendre. La solución de ésta 
ecuación se escribe como un serie de potencias, 



oo 

P(x) = ^ a^x^ 

j=0 



(9.25) 



donde el parámetro a no se conoce. Al derivarla da, 



dP{x) 
dx 



ax 



j=0 



3=0 



(9.26) 



y por tanto 



1 - x"^) — ^ax" 
dx 



ajx^ + x" jüjX^ ^ 



3=0 



3=0 

oo 



3=0 j=0 



(9.27) 



Derivando nuevamente y reemplazando en la ecuación generalizada de Le- 
gendre, 



d 

dx 



;i - x^) 



2,dP{x) 



dx 



m 



1 — x^ 



P{x) = a(a - l)x'^-^ J2 



3=0 



3=0 



+ l)a;" 



3=0 

oo 



3=0 



üjX^ — ax'^'^^ 'y^^jttjX^ 



3=0 



^ - (a + 2)a;"+^ ^ JCij^^'^ 

3=0 3=0 



— X 



a+2 



3=0 



-/(/ + !) + 



m 



1-x^ 



Y ^j^^ = O' (9-28) 

3=0 
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donde se verifican las ecuaciones, 



-/(/ + !) + 



X 



[2 



ax 



1 — X^ 

}oo 
^ OjX-' = O 
j=0 

- (a + 2)x"+^] 

oo 

^ jajX-?"^ = O 



3=0 



3=0 



(9.29) 



que al solucionarlas y al reemplazar los valores en la ecuación original, (ver 
|6], ¡9], [21] y [23]) se llega a una función recursiva para la variable aj de la 
forma, 



a,+2 



(a+j)(a + j + l) - ¿(¿ + 1) ^ 
(a+j + l)(a+j + 2) 



(9.30) 



lo que da la solución de las funciones Pj{x) en forma de polinomios, los cuales 
toman el nombre de polinomios de Legendre, que se compactan mediante la 
fórmula de Rodríguez como. 



P3Í^) 



1 é 

2H\ dx^ 



{x^-iy 



(9.31) 



9.3. ECUACIÓN LAPLACE EN CE 



87 



Algunos términos de la serie son: 



Po{x) 


= 1 




Pi{x) 


= X 




p (t^ 
-r2{X) 


- 2 ) 






= ^ {5x^ - 3x) 




Pi{x) 


= - (35x^ - 30x^ + 3) 
8 




P^x) 


= l (63x^ - 70x^ + 15a;) 
8 




Peix) 


= (231a:^ - 315a;^ + lOSx^ 
16 ^ 


-5) 


Plix) 


= — (429a;^ - 639a;^ + 315a;^ 
16 ^ 


-35a;) 



: (9.32) 

9.3.1. Aplicaciones en la Electrodinámica 

Una de las aplicaciones más usadas de los Polinomios de Legendre y las 
funciones de Rodrigues está en la expansión multipolar, que es necesaria en la 
teoría de radiación. Los polinomios ayudar a expandir funciones tales como, 

1 1 °° ^.'i 

1 TT = > = V -irrPiicos-f) (9.33) 

\x-x'\ ^r^ + r'^ -2rr'cos-f ^ r'^' 

para la restricción multipolar, r > r', la que es usada, por ejemplo, para 
calcular el potencial de una carga puntual, el cual usando el eje z como eje 
de simetría da, 

oo 

$(r, e)^Y^ [Ay + S,r-('+^)] P,(cos 7) (9.34) 
1=0 

donde Ai y Bi se hallan con las condiciones de frontera del problema. 

La expansión multipolar puede realizarse primero haciendo un expansión de 

Taylorde una función arbitraria T{r — R) alrededor del origen r = O de la 
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forma, 

^(r -R)^ r{R) + Y, ^í-^íí^) + ^ Z¡ ^-i^j^iÁR) + ■ ■ ■ (9-35) 



donde 



2 

i=x,y,z i,j=x,y,z 



3 / T-=0 

Si ^i(r — i2) satisface la ecuación de Laplace 



^^^•(^) - ( Q,.,. ) ■ (9-37) 



[V'T{r - R)] = 5] = O, (9.38) 



entonces la expansión se escribe tensorialmente de la forma, 

^ nr,:F,,{R) = \ Y, (3r,r,-¿,,r2) J-,,(ñ) (9.39) 



3 

i,j=x,y,z t,j=x,y,z 



si se considera J^{r — R) = |^rrR|' entonces J^{R) — ^, J^i{R) — 

;, defir 

polo y el cuadrupolo por 



^ij{R) = ^^2^2-p2^_^ etc, definiendo, respectivamente, el monopolo, el di 



N 



1=1 

N 

a=l 

N 



a=l 

: (9.40) 
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donde se obtiene la expansión multipolar de la función escalar, denominada 
potencial total $(-R) 

1 ^ 

1=1 

i=x,y,z 



+ ^ E Q^J{^R^RJ-kR')+■■■ (9.41^ 



6R^ 

i,j=x,y,z 

Si se aproxima a dos cargas, entonces el potencial dipolar eléctrico estará dado 
por, 

/ „N 1 P • R 

En coordenadas esféricas toma la forma (ver referencias ¡6] , ¡9] , [21] y [23] ) , 

^ oo / TV 

"^(^) = i^^E E i-ir'iriR)-Y.i^Kir.), (9.43) 

o 1=0 m=-l i=l 

donde se denomina un armónico sólido y está definido por. 



y RY^{ri) se denomina armónico sólido regular. Definiendo el momento mul- 
tipolar esférico como = ^^^iQíRTÍ'^í) para —l<m< —l, se halla el 
potencial como, 

^ 00 i 

HR) = ^ E E i-^ririR)QT 

o /=0 m=-l 



^ 00 

— E 



1=0 



4n 1 



Aneo ^ y 21 + 1 R^+^ 



J2 {-irYriRjQT (9.45) 



m=—l 



para R > rmax, y donde se observa que la expansión aparece como coeficientes 
en la expansión R~^ del potencial. 
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9.4. Clases de Gauge 

Algunos tipos de gauge para /i = 0, 1, 2, 3 y para j = 1, 2, 3 son el gauge 
de Lorenz, d^A^ — 0; el gauge de Coulomb ó de radiación, 'V .A — djA^ — 0; 
el gauge axial, A^ — 0; el gauge temporal ó Hamiltoniano, Aq — 0; el gauge 
de cono de luz, n^A'^ = 0; el gauge de Poincaré, XjAj = O y el gauge de 
Schwinger-Fock, X/j^A^ — O 



9.5. Campo de un Solenoide Largo y Delgado 



El campo de un solenoide largo y delgado se estudia desde el punto de 
vista del potencial vectorial de la forma. 



1 r J{^^,^, (9_4g) 

c J \x — x'\ 



en donde se hace la aproximación |a; — a;'| ^ = |a;| ^ + + . . . para obtener, 

^, , 1 fJ{x')^. , 1 [ {x.x')J{x') , 

Aix) = - ¡ -^d^x' + - / ^ , ^ ^ d^x' + ... (9.47) 

c J \x\ CJ \xY' 

que para un sistema se escribe de la siguiente manera: 

AAx) = ^ í Mx')d^x' + . í x'JAx')d^x' +... (9.48) 
c\x\ J c\x\'^ J 

Teorema. 

Si f{x') y g{x') son dos funciones bién comportadas de x'. Si J{x') está lo- 
calizada y tiene divergencia nula, entonces 

J ifJ . V'g + gJ . V7) d^x' = O (9.49) 

que se puede escribir como 

J [fV'.{gJ)]d'x'^0 (9.50) 
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integrando por partes mediante el uso de u = f, du = V'/á^x', dv = 'V' . 
(gj) d^x' , V = gJ y de U = g, dU = Vgd^x', dV = V . (gJ) d^x' y v = gj 
queda, 

J [fV . (gJ)] dV = - j{gJ. V'f)d'x' = O (9.51) 

pues fgJ evaluado en el universo es nulo. Con los valores f = 1 y g = x'^ 
entonces se tiene que V'/ = O y 'V'g = x[ y por tanto la integral queda, 

J . V'x'i) d^x' = j Ji{x')d^x' = O (9.52) 

Ahora si se hace la elección f = x[ y g = x'j entonces se tiene que V'/ = x[ 
y V'^f = x'j y por tanto la integral queda. 



J {fJ.V'g + gJ.V'f)d'x' 
= j {x'iJ .x'j + x'jJ .x[) d^ 

{x¡.Jj + x'j.J,;) d^x' (9.53) 



y por tanto. 



4.,) .V^ o (9.54) 

Se concluye, debido a la ec. (19.521) . que el primer término del potencial vec- 
torial (19.471) se anula y por tanto el término siguiente, usando la ec. (19.541) . 
se modifica de la siguiente manera. 



c\x 



3 



J J,{x')x'd^x' = . J x'Ji{x')d^x' 
Ji{x')d x' 



c\x 



3 Z^- 



3 



i)d^x' (9.55) 



debido a que x[Jj + x'^Ji = O, donde x-Jj = —x'jJi, entonces 2a;-Jj = 
— {x'jJi — x\Jj) y por tanto x[Jj = —\{x'jJi — x[Jj). Se nota que el resul- 
tado está relacionado con el delta de Levi-Civita, por lo que el término se 
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modifica a 



^ ^ ] ^-^i / (•^j'^i XjJi)d X 

ijk 

^— --a; X I (x' X J(x')(fx' 



2c| 



1 m X a; ^ , , , , 

a; X m = — — --— = A{x) (9.56) 



donde m es el momento magnético y /x(x) = es la magnetización del 

medio. Entonces la densidad de flujo magnético está deflnida por, 



„ , m X X , , 

B = V X A = V X , . 9.57 

\xr 



Usando 



[V X (A X B)]. = EijkVjSklmAlBm 
^ ij klm^ j Al B fn 

= {SaSjm - SimSji) VjAiBm (9.58) 

donde usando las propiedades del los delta de Kronecker se obtiene: VjAíBj — 
V jAjBi = (VjBjAi + Bj'VjAi) — (VjAjBi + AjVjBi), dando como resultado 
la ecuación vectorial, 

V x{Ax B) = {V . B)A + {B . V)A - (V . A)B - {A . V)B (9.59) 

que al aplicarla en la ecuación (19.571) da el siguiente resultado: 

„ 3x (x . m) — m , ^ 

B = ^— ^ 9.60 

\xr 



Y así, el potencial vectorial es (en el sistema CGS-Gauss), 
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Si m está alineado en la dirección del eje z, el potencial en el punto P{r, 9, 0) 
es, 

. , , 1 m sin ^ - , , 

A{x) = —(j) (9.62) 



y por tanto la densidad de flujo magnético se puede calcular mediante, 



c 



1 d f m sin 9 



r sin 9 89 
dando finalmente. 



sm I 



1 d f msm9 



r dr 



r e 



— {2cos9f + sm9d 



(9.63) 



(9.64) 



que es similar al campo eléctrico dipolar, con momento dipolar eléctrico p: 

E ^^{2cos9r + sin9e) , (9.65) 

pero a pesar de dicha afinidad, el potencial escalar no se usa para representar 
el campo magnético dipolar. Para un solenoide muy delgado con momento 
magnético por unidad de longitud M.{x') a lo largo de su longitud, el po- 
tencial vectorial en un punto distinto de su eje de simetría longitudinal es. 



A(.) = - 



1 f M{x') x{x- x') 



'13 



M 

c 



\x — X 



dL{x') 



(9.66) 



donde r' — -\- z"^ — 2 zr' eos 9 (Ley de cosenos) y sin^' = ^sm.9 (Ley de 
senos), y por tanto el potencial vectorial en un punto diferente a lo eje de 
simetría queda 




M{x') x{x- x') 



\X — X 



,/|3 



dL{x') 



sin 9' 



L (r^ + z"^ — 2z'r eos 9) 2 
sin 9' 

1+ (^)^ -2^cos^ 




(9.67) 
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pero r~'^z"^ — 2r'^ cos9 + l = {r~^z' — cos6) — cos^ 6 + 1 = {r~^ z' — eos 6) + 
sin^ 9. Haciendo en cambio de variable u — r~^z'— eos 9, entonces du — r~^dz', 
y por tanto, 



A{x) 



M sin 9 



rdti 



-(^+cose) {u'^ + sm'^9) 



4> 



(9.68) 



realizando el reemplazo trigonométrico u = sin 6* seca, se resuelve la integral 
resultando. 



A(x) 



M 



c rsin9 



eos 9 + 



^ + eos 61 



+ eos 9)^ + sm'^9 







(9.69) 



si se considera la longitud L del solenoide, el ángulo 9 formado entre el eje 2; y 
un vector r que va desde el punto superior del solenoide al punto de medición 
P, y el ángulo 6*2 formado entre el eje 2; y un vector r2 que va desde el punto 
inferior del solenoide al punto de medición P, se observa que la ecuación se 
reduce a, 

A(x) = — (eos 92 - eos 9) 4> (9.70) 

c r sin é* 

dando el valor del potencial vectorial en un punto diferente a lo eje de si- 
metría. Por tanto la densidad de flujo magnético será. 



B{x) 



Vr X A{x) + X A{x) 

Id 1 d - 

—— — (A^sm9)r - - — {A^r)(f> 

rsm.9o9 r or 



+ 



1 



d 



r2 sin 92 d92 



{A^ sin 92) f 



resultando 



B{x) = —r 2^2- 



(9.71) 



(9.72) 



El resultado anterior es una de las razones para la existencia del monopolo 
magnético, pues tien la forma de dos polos magnéticos puntuales con cargas 
magnéticas positiva y negativa, donde la magnitud de éstas es 
resultando (para los puntos por fuera del eje de simetría). 



M 

c 
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con un potencial vectorial magnético dado por, 

A(x) = (eos 02 - eos 6) 0, (9.74) 

r sin 9 

donde se puede ver una correspondencia con la parte eléctrica tanto del po- 
tencial escalar (f){x) = - — — , como de la intensidad de campo eléctrico 

En el caso de que las ecuaciones fl9.70p y (19.721) sigan la restricción 6*2 = O 
para L ^ oo o sea para r2 — oo (con el fin de hacer una aproximación a un 
extremo puntual), quedan como 



7->/ \ Qm ^ 

B{x) = —r 

A(x) = -SlIL^(l-cose)(f> 
r sm O 



que no están definidas sobre el eje z negativo, debido a que el solenoide 
delgado o magneto está sobre dicho eje, que representa la distribución de 
corriente fuente. Pero un monopolo real tendría un campo esférico isotrópico 
B{x) = que no lo mismo que B{x) = ^f, para 6 ^ n. Se observa que 
usar un solenoide semi-infinito para representar a un monopolo magnético 
no es correcto, aunque la línea de singularidad tiene un significado muy claro 
para el solemoide. Para el monopolo ha sido creado el potencial vectorial. 



r sin 9 r \2 

con el objetivo de reemplazar al solenoide semi-infinito por un monopolo 
puntual, a pesar de sus inconsistencias físicas. 



9.6. Principio de Mínima Acción 

El principio de mínima acción o principio de Hamilton se basa en que para 
un sistema que evoluciona en el tiempo, la diferencia entre su energía cinética 
y su potencial, para la trayectoria mas realista, es mínima. La diferencia entre 
estas energías se denomina la función Lagrangiana, £ = T — ?7, y la acción, 
S, es la evolución temporal de ésta Lagrangiana, 

/"Í2 

S= C{q,q,t)dt (9.76) 
Jti 
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donde g, q son las posiciones y velocidades generalizadas del sistema. Según 
el teorema fundamental del cálculo variacional, la variación de la acción debe 
ser mínima y a primer orden en la trayectoria y se debe anular en los puntos 
extremos, lo que genera una integral extremal. Por tanto, 

(9.77) 



SS^S Í £{q,q,t)dt 
Jti 



lo que genera variaciones de las coordenadas generalizadas q{t) — * q{t) +6q{t). 
Adicionalmente, debido a la integral extremal se tiene que Sq{ti) — Sq{t2) — 
O, entonces. 



Í2 



SS^S C{q + Sq,q + Sq,t)dt-S / C{q,q,t)dt 



ti 



Í2 



(9.78) 



como el principio variacional exige que la primera variación sea cero entonces. 



ss 



Í2 



dC ^ 9£ _ , 
—dq + —5q ]dt^O 
oq oq 



al integrar el segundo término por partes se llega a. 



'—Sq 
dq 


Í2 

+ 

íl 


rt2 


'dC 


d 




Al 


dq 


dt 





5qdt = O 



(9.79) 



(9.80) 



íl 



pero debido a la condición de extremal se tiene que: 

debido a que las variaciones de la Lagrangiana son independientes de Sq que 
es distinta a cero, entonces se verifica. 



Í2 



= 0. También 



d /dC\ _dC _^ 
dt \dq J dq 



(9.81) 



Si hay varios grados de libertad, entonces deben existir diferentes acciones, 
y por tanto variaciones de éstas para las s funciones diferentes qi{t), lo que 
lleva a 



d_ íd£\ _ ^ _Q 



dt \ dqi J dq 



(9.82) 



para i — 1, 2, 3, . . . , s, denominadas ecuaciones de Lagrange. De éstas se 
observa que si el potencial no es dependiente de la velocidad. 



d [ d[T{qi) - Ujqj 

dqi 

dqi 



dt \ 

d 



d[T{q,) Ujqj 
dqi 



O 



'dT\ 
dt \dqij 



(9.83) 
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si pi = 1^, es el momento canónicamente conjugado, entonces: 

dpi dU 

Í = ^- = -aí 

que verifica la dinámica de Newton. 

Además si se calcula la variación temporal de la Lagrangiana (usando suma 
sobre índices contraidos), entonces, 

dC dC doi dC dái dC ,^ 
= 11 _| ±L _| (9.85) 

dt dqi dt dqi dt dt 

aplicando las ecuaciones de Lagrange (19.821) 

dC d í dC\ , dC dqi dC 
dt dt \ dqi ) ^* dqi dt dt 
d ( dC\ dC 



y por tanto se verifica que 

dh dC 
~ctt ~ ~~dt 



(9.87) 



donde h = qi^ — C. Ésta función h{qi,qi]t) se denomina función energía y 

si la Lagrangiana no depende del tiempo, entonces se conserva. Esta trans- 
formación es un transformación de Legendre, y si se incluye el momento 
conjugado Pi = entonces el diferencial de la Lagrangiana es 

dC 

dC = pidqi + pidqi + —dt (9.88) 

ot 

y se define la función Hamiltoniana mediante una transformación de Legen- 
dre, de la forma 

n{q,p,t) = qiPi-Ciq,q,t) (9.89) 
con su diferencial igual a 

dn = dqiPi + ¿lidpi - ^ 

dn, dn, dn,^ 

= -g^dqi + —dpi + —dt (9.90) 

donde la Hamiltoniana no depende de las velocidades generalizadas y por 
tanto la transformación de Legendre produce un cambio del espacio de las 
configuraciones {qi,qi) al espacio de las fases {qi,Pi) 
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